Theory of minimum spanning trees II: exact graphical methods and perturbation 

expansion at the percolation threshold 
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Continuing the program begun by the authors in a previous paper, we develop an exact low- 
density expansion for the random minimum spanning tree (MST) on a finite graph, and use it to 
develop a continuum perturbation expansion for the MST on critical percolation clusters in space 
dimension d. The perturbation expansion is proved to be renormalizable in d = 6 dimensions. We 
consider the fractal dimension Dp of paths on the latter MST; our previous results lead us to predict 
that Dp — 2 for d > dc — 6. Using a renormalization-group approach, we confirm the result for 
d > 6, and calculate Dp to first order in e = 6 — d for d < 6 using the connection with critical 
percolation, with the result Dp = 2 ~ e/7 + 0{e'^). 



I. INTRODUCTION 

In this paper, we continue our treatment of the statis- 
tics of random minimum spanning trees, which was be- 
gun in Ref. [l| (to be referred to as I). We first recall the 
definitions (see any of Refs. 0, H, S IB ) ■ '^^ consider an 
undirected, connected graph G with vertex set V, edge 
set E and a real-valued cost ie assigned to each edge 
e & E. A spanning tree is then defined as a subset of the 
edges of G that connects all the vertices and contains no 
cycles: in other words, it is a tree and it spans V. Such a 
tree must exist because the graph is assumed connected. 
A minimum spanning tree T is a spanning tree such that 
the sum of the costs of its edges. 



(I.l) 



eGT 



is minimized over the set of all spanning trees on G. (If 
the costs ie are strictly positive, then any spanning sub- 
set of the edges that has minimum cost is automatically 
a tree.) In this paper, we again consider the random 
problem, in which the edge costs are assumed to be in- 
dependent and identically-distributed (iid) random vari- 
ables, with a continuous probability distribution for the 
cost of each edge, and we are interested in the statis- 
tical geometry of the tree. This random model will be 
referred to simply as the MST. We also consider a gen- 
eralization, introduced in Ref. T|, in which the costs are 
iid, the probability distribution for the cost of any edge is 
uniform on the interval [0, 1], and one finds the minimum 
spanning forest (a forest is a collection of vertex-disjoint 
trees) on the (not necessarily connected) graph formed 
by the subset of edges that have cost less than p G [0, 1]; 
this object is called MSF(p). It is connected with the use 
of Kruskal's greedy algorithm for the MST [1, i, i, i, Q . 
In Kruskal's algorithm, the edges are tested one by one, 
in order of increasing cost. Each edge is "accepted" as 
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belonging to the MST unless doing so would form a cycle 
in combination with edges already accepted. If this algo- 
rithm (or "Kruskal process") is stopped when all edges 
with cost < p have been tested, one obtains the MSF(p). 
MSF(p) is closely related to bond percolation which 
is obtained if the process is modified so that it simply 
accepts every edge, in which case each edge is indepen- 
dently accepted (occupied) with probability p. This re- 
lationship plays a central role in the work to be reported 
here. 

The background and motivation for studying this prob- 
lem were discussed in the previous paper. Of particu- 
lar interest there was the strongly-disordered spin-glass 
model of Newman and Stein (NS) Our results 

strongly suggested that when the graph is a lattice in 
Euclidean system of dimension d, the critical dimension 
for that model and for the MST is 6. This is the dimen- 
sion above which there are many large connected compo- 
nents of the MST that are visible within a window of size 
W; we found that the number is of order W^^^ . The 
strongly disordered spin glass maps onto the MST, and 
the logarithm of the number of ground states of the spin 
glass that can be distinguished within a window is the 
same as the number of large connected components that 
are visible. For the MSF(p), the properties such as the 
number of large connected components are expected to 
be similar (scale with the same exponents) for all p > pc, 
where Pc is the threshold for bond percolation. These re- 
sults were strongly motivated by the solution of MSF(p) 
on the Bethe lattice (Cayley tree) 

Other problems related to the MST that have been 
studied in the physics (and related) literature include the 
fractal dimension of the paths on the MST in Euclidean 
space, and the relation with problems of optimal paths 
and transport in random media 0, [loj . In this paper we 
begin to address some of these issues. On a finite graph, 
there is a unique self-avoiding path on the MST between 
any two vertices. Eventually, we aim to study the frac- 
tal dimension of the paths on the MST. In general, 
such a dimension can be defined by a box-counting tech- 
nique, that is counting the number of steps on the path 
inside a window. One would hope also to obtain the frac- 
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tal dimension of the path from the expectation value of 
the total number of steps on the path, if this is of order 
i?^p for two vertices separated by a large Euclidean dis- 
tance R. [In the Kruskal process or MSF(p), this would 
again be expected to be the same (on large length scales) 
for any p > pc, for two vertices that are connected on 
MSF(p).] Several earher works [ll|, [H [H, El have pro- 
vided numerical estimates of D'^ for the MST at p = 1 
in various spatial dimensions. In rf = 2 dimensions, the 
value is around 1.23. 

However, in the limit of an infinite system in d > 6, 
there are many connected components of the MST, and 
for a large finite system (say, a hypercube of side L) this 
means that for many pairs of vertices, the path connect- 
ing them makes a large excursion (of order the system 
size). To see this, first note that if one makes the system 
larger by adding additional edges at the boundary, cho- 
sen from the same probability distribution, keeping those 
in the interior the same, then one can study the behavior 
of the MST as L increases. Locally, it will converge to 
a definite limit that has the number of components as 
described above. This implies that the path connecting 
two fixed vertices will be deformed until it goes off to 
infinity as L —^ oo, and so have infinite length, unless 
the two vertices are on the same component in the limit. 
For finite L, it will make an excursion typically of linear 
size of order L. These considerations, together with the 
probability of order l/R^^^ that the two given vertices 
are on the same connected component in the limit, sug- 
gest heuristically that the expected length of the path 
scales as 

0(i?^p /R'^-^) + ©((!- l/R'^-^)L^'^^ (1.2) 

where D'^ is the fractal dimension of the path on the 
MST, for which the Bethe lattice results in Ref. [l| sug- 
gest = 2 for d > 6. For L > i?, the last term 0{L'^) 
dominates, and again the result is similar for MSF(p) 
for all p > pc- But if we could restrict attention to (or 
condition on) paths not going to infinity, the conditional 
expectation for the path length would scale as R^p . For 
d < 6, this problem does not arise, and we expect that, at 
least for d < 6, where there is presumably only one con- 
nected component in the infinite system limit, the length 
of the path will indeed scale as R^p , independent of L 
for large L. (In the numerical work cited above, cases in 
which d > 6 were apparently not considered.) 

For this and for other related technical reasons, we will 
restrict our calculation of exponents in this paper to the 
Kruskal process with p < Pc- In this case the probability 
that two points are connected by either the percolation 
or the Kruskal processes decays with increasing distance 
R. For p = Pc, where scaling will again apply, we will 
denote the fractal dimension of the path on the MSF(pc) 
by Dp. The naive scaling for the total length of a path 
R^p still works in this case. 

There are claims (the "superhighways" argument) 
^Ol [ist that in the Kruskal process (on various fami- 



lies of graphs) the properties of the paths are mainly 
determined by the percolation threshold. Numerically, 
the fraction of steps on a path on the MST that are al- 
ready present at p = Pc, averaged over all pairs of end 
vertices for the path, goes to a constant [Tol . [T5| . This 
suggests that D'^ = Dp. For the Bcthe lattice and for 
Euclidean systems with d > 6 we can argue that both di- 
mensions equal 2 so equality holds, however, it is less 
clear whether it holds for d < 6. Some support for it can 
be obtained from scaling arguments in percolation. Con- 
sider a hypercube of side L, and for bond percolation ask 
for the probability that two given opposite faces are con- 
nected when the occupation probability for each edge is p. 
The percolation threshold is the value of p above which 
the connection occurs with probability one as the size 
L — > oo. For L finite, connection occurs at p = Pc with 
low probability, but it occurs with probability approach- 
ing 1 at p — pc of order l/L^^'^ , where ly is the correlation 
length exponent in percolation J|. At this value of p, the 
correlation length ^ is of order {pc — p)^" — 0{L), and 
scaling properties at scales less than ^ should be the same 
as those at threshold p = Pc. Hence we expect the frac- 
tal dimension Dp of the path on the MST connecting the 
two faces (which is one of the paths on the percolation 
cluster that do so when the faces are first connected as 
p increases) to be Dp. However, this path does not have 
the boundary conditions we wanted, as we only asked for 
the connection of the two faces, not of two given vertices 
in the interior of a system. When the separation R of 
the vertices is large, we may surround each by a nested 
sequence of concentric spheres of radii say R/3 for 
j = 0, 1, ... and ask the same question for each annulus 
bounded by two of these spheres. Then the relevant p 
will be different for each sphere (higher p is required to 
make the connections to the vertices on smaller scales), 
but also the scaling holds for each one. In addition, the 
clusters that connect each pair of spheres must also be- 
come connected together to form a single cluster and an 
MSF path. Similar arguments apply to all of these. This 
does suggest, heuristically, that Dp = Dp. 

On examining this argument, a key part of it can be 
seen to be the idea that there is a unique candidate su- 
perhighway (critical percolation cluster) that is used to 
make connections over large distances on each scale, and 
so it is clear which ones must be connected by higher- 
cost "roads". This is a property of critical percolation 
clusters that holds for d < 6, but not for d > 6, where 
the number of large clusters visible in a window of size 
W is W'^~'' 0,[il]- This behavior may itself underhe the 
result [l| that the number of connected components of 
the MST (or of MSF(p) for p > pc) has this same form 
W^-^ for d > 6 (but order 1 for d < 6). 

In this paper we do not assume the equality of Dp 
and Dp, but will study Dp. We first construct in sec- 
tion |TT] an exact series expansion for any finite graph that 
gives the probability that the path on MSF(p) from ver- 
tex X to vertex y passes through vertex z, by analogy 
with expansions in percolation. The expansion takes the 
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form of a weighted sum of subgraphs. This expansion 
may be of general interest within various approaches to 
MSTs not considered further in this paper, such as high- 
temperature series. 

In section IIII[ we then turn this expansion into an 
asymptotic or perturbation expansion in the continuum 
(with cutoff) by neglecting excluded volume requirements 
that were present in the exact expansion, and taking the 
graph to be the infinite lattice. (From here on, our results 
are not fully rigorous mathematically, though they will 
satisfy most theoretical physicists.) In this way, we ob- 
tain a Feynman diagram expansion. It is related to that 
for percolation, but with modified Feynman rules. In or- 
der to avoid dealing with the appearance of the order 
parameter for percolation (non-zero probability of con- 
nection to infinity), from this point forward we consider 
only p < Pc- 

The perturbation expansion contains ultraviolet 
(short-distance) divergences as the wavevector cutoff goes 
to infinity (or as the lattice spacing goes to zero) with the 
separation of the vertices held fixed. We prove that the 
expansion is renormalizable in the field-theoretical sense 
to all orders in perturbation theory. We then use stan- 
dard techniques to formulate a renormalization group 
(RG) approach which gives the scaling behavior of the 
correlation functions (probabilities) already mentioned. 
The exponents or fractal dimensions are then calculated 
for d < 6 as an asymptotic series in powers of e — 6 — d, 
with the result 

Dp = 2-e/7 + 0{e^) (1.3) 

as e — > 0. For d > 6, the analysis of the cut-off expansion 
confirms that Dp = 2 to all orders in perturbation the- 
ory, and the path behaves as a Brownian random walk 
on large length scales. Further discussion is contained 
in the Conclusion. To improve readability of the paper, 
many detailed derivations have been relegated to the Ap- 
pendices. 

Our calculations can in principle be extended to other 
exponents, such as those defined in (l7j . or carried to 
higher orders in e. They can also be extended to in- 
clude statistical properties that involve the cost of the 
MST. There do not appear to be any scaling relations 
that relate the geometric exponents for MSTs to those 
for percolation, unlike those found for the costs in p^ . 
even though the critical dimension dc = 6 is the same. 



II. LOW-DENSITY EXPANSION ON A FINITE 
GRAPH 

In this section we describe how to set up an exact low- 
density expansion which enables one to calculate con- 
nectedness functions for the (random) MSF on any finite 
graph G. The expansion takes the form of a low-density 
expansion, similar to a high-temperature expansion fa- 
miliar from statistical mechanics, with small p playing 



the role of high temperature, and is modeled on a corre- 
sponding expansion for bond percolation. Although from 
the point of view of this paper the primary utility of the 
result is to provide a basis for the RG calculations of the 
following section (which determine the fractal dimension 
of paths on the MSF), the expansion is of interest in its 
own right, and might form the starting point for further 
mathematically-rigorous calculations, using e.g., lace ex- 
pansion methods [l9|. As the construction is somewhat 
involved, we concentrate here on describing the structure 
of the expansion, and relegate most details of its deriva- 
tion to the Appendices. 



A. Graphical expansions for percolation 

In view of the correspondence between bond perco- 
lation and Kruskal's algorithm discussed in I, and the 
relation of the MSF(p) process with percolation, we will 
first briefly review the low-density expansion for bond 
percolation. Afterwards, we set up a corresponding ex- 
pansion for MSF(p) by generalizing the arguments using 
Kruskal's algorithm. 

Bond percolation on a finite-dimensional lattice is con- 
ventionally treated as the Q — > 1 limit of the low-density 
expansion of the Q-state Potts model 0, [H [H, HI . The 
partition function of the Potts model, which is a polyno- 
mial in Q, is also known as the Tutte polynomial |24| . 
and is a generating function for Q-colorings of the ver- 
tices of the graph G with weights that depend on whether 
adjacent vertices are given the same or different colors. 
Although the Q ^ 1 limit of this function can be taken, it 
lacks a mathematical definition in terms of state variables 
(colors), so in order to establish a correspondence with 
percolation we instead use a method originally due to 
Essam [25l . [26j . This has the advantage of being phrased 
explicitly in terms of geometric quantities. 

The basic object of interest is the two-point connected- 
ness function Cx,y(p), the probability that two vertices x 
and y on the lattice are connected by a percolation clus- 
ter (connected component) when the probability that an 
edge is occupied is p. In Appendix lA H we review the 
graphical expansion for this probability, which has the 
form 

Cx,y(p)= d{G)Fr[G<pl (II.I) 

where the sum is over all graphs G in the set Gx.y, de- 
fined as graphs on (i.e. subsets of) the lattice containing 
the endpoints x, y. The factor of Pr[G' < p] is the prob- 
ability that all the edges of G have cost less than p (i.e. 
are occupied in the percolation process), which in the 
iid model (Bernouilli model of bond percolation) is sim- 
ply (where E{G) is the set of edges of G). The 
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function d{G) can be defined as 

d{G)= J2 connects x,y]; 

E'CE{G) 

(II.2) 

where the sum is over subsets E' of the edge set E{G). 
We use the notation I[X] for the indicator function on 
events X, which takes the value 1 when X is true and 
when X is false. d{G) = if a; and y are not connected 
byG. 

The preceding graphical expansion may be general- 
ized straightforwardly to n-point connectedness func- 
tions: now we must sum over Gxi,- -,x„, the set of sub- 
graphs of the lattice that contain all n vertices, and in 
Appendi x I A 21 we show that the correct generalization of 

(iir^ is [23 

E'CE{G) 

X I [E' connects xi, • • • x„] . (11.3) 

The d-weights defined above have several properties 
worth remarking on. First, in Appendix lA 3[ we prove 
that the expansions pi.2p , (jlLSp assign the same weights 
to the same diagrams as does the Q ^ 1 limit of the low- 
density expansion for the Q-state Potts model, so the 
above expressions are entirely equivalent to the conven- 
tional description of bond percolation. 

Second, the d-weight clearly only depends on the con- 
nectivity of the set E{G) of the edges of G, not on the 
geometry of G, and further as shown in Appendices lA 11 
lA 21 the d-weight is invariant under replacing edges of 
G with chains of edges (i.e., inserting vertices of degree 
two). Thus d{G) is a topological invariant of graphs 
G with two marked vertices. This means that we can 
consider the expansion in terms of topological graphs G, 
which are simply graphs (without any embedding in the 
lattice), with two distinct vertices labeled x, y (and the 
others unlabeled) , and which may be assumed to contain 
no vertices of degree two other than possibly x, y. Using 
this property, we may rewrite the expansion pi.ip as 

C^Ap)= E 4m ^ Pr[G<p], (II.4) 

with analogous expressions for the n-point functions. 
Here C/x,y is the set of all (topological equivalence classes 
of) topological graphs with two distinct labeled vertices 
X, y, and the inner sum is over all possible embeddings 
X : Q —> G which map the edges of Q into self-avoiding 
chains of edges on the lattice, producing the set of lat- 
tice graphs summed over in (jll.ip . Note that these chains 
must not only be self-avoiding walks, but also must avoid 
intersection with chains arising from different edges of G- 
If the topological graph G has any non-trivial automor- 
phisms (leaving the root points x, y fixed), then there is 
more than one way to produce the same embedded graph 
G. Consequently, we must divide by the number A{G), 



which is the number of elements in the automorphism 
group of G- This number is often called a "symmetry 
factor" , especially in the context of Feynman diagrams, 
and will play such a role later. 

This topological property of the d-weights is crucial 
for extending the lattice expansion (|II.4p to a continuum 
theory, a point to which we will return in Section Fill Al It 
also simplifies lattice calculations, since it greatly reduces 
the number of different graphs for which d(G) must be 
calculated. The d-weights have further properties that, 
when the graph G has a connected subgraph, allow them 
to be factorized into pieces given by the d-weight of the 
subgraph and that of the "quotient" graph in which the 
subgraph is replaced by a single vertex. These are dis- 
cussed in the Appendices. 

There is a further function will be useful in making 
comparisons with the MSF theory. This is the derivative 
of Gx.y(p) with respect to the value pe oi p on a particu- 
lar edge e (the generalization of the formulas to cases in 
which the occupation probabilities such as Pe for edges 
differ should be obvious). The probability that x and y 
are connected by a cluster at p changes with pe only if x 
and y are not connected when pe = 0, and are connected 
when Pe = p. This implies that for p^ = p, any path from 
X to y on the cluster must traverse e. An edge with this 
property is called a singly-connected edge. So we define 

G^,y(p)= Tdp ACx.y(p,P,), (II.5) 

Jq Ope 

which is the probability that at parameter p, x and y 
are connected, and e is a singly-connected edge on the 
same cluster. Note that C^y{p) is not the same as the 
3-point connectedness function that was defined above, 
and its lattice expansion (which may be obtained directly 
from the preceding definition) still contains the same d 
weights as for Gx,y(p), and in fact is given by the same 
expansion (jll.4p . but with the additional condition that 
e be an edge of the embedded graph G (note that for 
the graphs G in the expansion, e does not have to have 
the singly-connected property). We may choose to view 
the topological graphs as having the inverse image of e 
as a marked edge (either of the ends of which may be 
degree- two vertices), so that the embeddings A map it to 
the single edge e; in this case the relevant automorphisms 
of G must fix this edge as well as x, y, and we denote 
the number of these by A'{G)- Clearly A'{G) < A{G) 
(one group of automorphisms is a subgroup of the other). 
These different ways of writing the function are equiva- 
lent. 

Although we have formulated the expansion here in 
terms of the infinite lattice, it proceeds in exactly the 
same way if the lattice is replaced by any finite graph 
G. Indeed it is best viewed as derived from some finite 
graph such as a portion of the lattice, followed by an in- 
finite volume limit. The sums over all embedded graphs 
make sense for a finite graph because only a finite number 
of terms contribute. For the infinite lattice, the sum con- 
verges for p < pc (like a high-temperature expansion). 
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but not for p > pc- In the latter case, it needs to be 
re-summed, but we will not discuss this here. 

In the next two subsections, we proceed to develop 
a graphical expansion for MSF connectedness functions, 
analogous to (|II.4p in that it takes the form of a weighted 
sum over topological graphs and their lattice embeddings. 
This will be the basis for the continuum theory analyzed 
in section Ulll 



B. A low-density expansion for MSF paths 

We will study the random geometry of paths on the 
MSF by introducing a connectedness function y (p) , 
which is the probability that two vertices x, y on the MSF 
are connected and the connecting path passes through a 
third vertex z. In the following, a vertex z satisfying 
this definition will be called a "MSF path vertex" (with 
respect to two other vertices x, y). (We could equally 
well define a similar function in terms of the probabil- 
ity that the MSF path passes through an edge e instead 
of a vertex z, which makes the analogy with the per- 
colation function C^y{p) closer; in the continuum for- 
mulation developed afterwards, there is no difference be- 
tween the vertex and edge cases for MSF.) We obtain 
a diagrammatic expansion for these connectedness func- 
tions by relating the Kruskal process defining MSF(p) to 
the expansion for bond percolation obtained in the previ- 
ous section. The difference between the percolation and 
Kruskal processes is that, in the latter, as p increases an 
edge with cost £ = p is not accepted if, together with 
edges already accepted hy i < p, it forms a cycle. 

In Appendices [Bl and O we develop the tools needed to 
obtain a diagrammatic expansion for the MSF connected- 
ness functions involving the MSF path vertex. Appendix 
IB] contains a careful discussion of properties of the MSF 
paths that are used, while Appendix [Cl uses the method 
of inclusion and exclusion to obtain the expansion itself. 
Here we will begin by defining notation. The ordering 
of the costs of the edges on a subgraph G of the given 
graph G plays an important role. We know that in fact 
the MST on a graph depends only on the ordering of the 
costs P, • We define an ordering on the set of edge 
costs of G as a permutation vr G S^e\ on the set of \E\ 

elements. We index the edges of G arbitrarily, and take 
the ordering of their costs to be defined by n via 

< 4o) -i < 3, (II-6) 

for the edges indexed i, j (writing £i for fgj. For subsets 
of E, such as E' or E{G), we define the induced order- 
ing by restriction, and denote it by tte' or tte(g)- Thus 
tte' is a permutation of the subset E' . We can obtain 
the probability for each possibility ordering from the iid 
probability distributions on the edge costs in the obvious 
way, and clearly Pr[7r£;'] = in particular for the 

special case E' = E. 



With these definitions, the result we obtain in Ap- 
pendix [U] for the path vertex connectedness function can 
be written 

X E dMSF(G|7r£;(G))Pr[7rB(G)], (II. 7) 

"■eco 6S'|E(G)| 

which should be compared with the corresponding result 
(|II.ip for percolation, or its analog for G^y{p). Here 
Gx,y;z is the set of graphs G that contain the root points 
X and y, and a (self-avoiding) path from x to y passes 
through z. The diagrammatic weight dMSF(G|7r) implic- 
itly depends on x, y, z, and is ordering-dependent. It 
can be defined as 

dMSF(G|7r) — 

J2 (-l)l^^(G)hlB'li[i;' connects x,y] 

E'CE{G) 

X I[7mst(G_b/ Ittb') passes through z], (II. 8) 

which should be compared with (jll.2p . Here Ge' is the 
subgraph of G with vertex set V and edges £", and 
7mst(G£;' Itt^;/) is the path connecting x, y on the MST 
on the graph Ge' with costs on E' induced by restriction 
from G. This MST can be assumed to exist because Ge' 
can be assumed to be connected, in view of the indicator 
function 1[E' connects x, y]. 

Next we wish to express this connectedness function 
as a sum over topological graphs G- We consider a given 
lattice graph G which is the image of some topological 
graph Q under some embedding A. For clarity in what 
follows, we denote elements of the edge sets of G, Q by 
different symbols: we have e S E{G) and e £ E{Q). We 
find it useful to extend our notation and let A(e) C E{G) 
denote the chain (self-avoiding path on the lattice) of 
edges on the lattice that the topological edge e is mapped 
to. 

For both percolation and MSTs, the only edge cost 
information relevant for the connectedness functions is 
the cost of the most expensive edge on A(e), which we 
denote by 

= max 4. (II.9) 

e6A(e) 

Then in bond percolation at occupation probability p, to 
determine whether A(e) connects its endpoints, we only 
need to check whether < p. Likewise, A(e) is a subset 
of the MST on G only if is less than the maximal edge 
cost encountered on all other paths on Q connecting the 
same vertices. We may think of as the cost induced on 
the edge e of by the embedding A. A corresponding or- 
der i^'E(g) is induced on these costs. Given an embedding 
of Q, both Pr[i(:] and Pr[7r^^g.^] depend on the embed- 
ding. The probabilities can be readily calculated, as we 
will see shortly. 
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Because dMSF(G|7r) is computed in terms of connected- 
ness properties (whether or not the MST path between 
the root points of the diluted graph Ge' goes through 
the MSF path vertex at z), it may be computed using 
only the relative ordering i^'^^g) obtained from the {^e}- 
Therefore we have 



dMS,F{G\TrE{G)) = dMsAQWEig)) 



(11.10) 



(and depends on the marked vertices of Q that are the 
inverse images of x, y, and z under A). In the expansion 
(III.7p . only the ordering and its probability is required, 
and it is possible to fix an ordering before choosing an 
embedding (and finally summing over both). Hence it 
may be written as 

C^x,y(p)= E E {dMSF{QWE(g))l^{9)) 

X E Pr[7r^(g) A(G<p)]. (11.11) 

Here again A'{G) is the relevant symmetry factor as de- 
fined in section III Al (strictly, it is defined here for auto- 
morphisms fixing the inverse image of the vertex z rather 
than e; the cases relevant to the continuum expansion 
later are those in which the vertex z has degree two, 
and should be compared with those in which, when e is 
shrunk to a single vertex, then it has degree two, and for 
these the numbers are the same). Eq. pi.ll[) is the main 
result of this Section, and should be compared with the 
percolation result pi.4p . It remains to find an expression 
for Pr[n'^f^g^ A (G < p)]. This is done in the following 
Subsection. 



This is clear from the iid behavior of the edge costs: the 
most costly edge could be any of the N2 edges among 
the total number Ni + N2- [The result is more general 
than the model we have been using, in which costs are 
uniformly distributed on [0,1]. For a general iid distri- 
bution of costs, the probability that Li < L2 < Iq is 
N^p^i+N^ /(jVi iV2), where p is the probabihty that a 
given edge is less than Iq- Similar statements apply to 
the following generalization.] 

The result generalizes to any embedded graph G, with 
edges E{Q) = {ei, . . . , e„} which have lengths A^i, . . . , Nn 
and maximum costs Li, . . . , L„. We let the ordering be 



'^Eig) = ^^^^ ^^^^ -^TT'Cl) < i7r'(2) < 

The generalization of pi.l3p is 



in)- 



Pr[7r' A (G < p)] = 



(n-l) ■ ■ ■ 



(11.14) 



which can be evaluated as 



Pr[7r' A (G < p)] = p 



1^1 Ej=i N^'U) 



= Pr[G<p]n 



(11.15) 



=1 Ej=i N^'U) 



C. Expression for Pr[7r^^gj A (G < p)] 

To put the quantity Pr[7r^^g,^ A (G < p)] in a more 
tractable form, we return to basic considerations. As 
the simplest example, we take the case where G consists 
of two root vertices connected in parallel by two edges 
ei, £2. We consider an embedding A where A(ei), A(e2) are 
chains of Ni, N2 lattice edges, the most expensive edges 
of which have costs Li,L2, respectively. We can define 
and evaluate the probability for Li < l\ and L2 < l2'- 



Pr[(ii < £1) A (L2 < ^2)] 



— K.^ K.2 



(11.12) 



The probability density for (Li, L2) is found by differen- 
tiating on both variables, and then the probability that 
Li < L2 < p is found by integration: 

Pr[Li <L2<p] 

di2 I dh——PG{ilJ.2) 



dti di2 
N2 



N1+N2 



M1+N2 



(H.13) 



Although the preceding result is completely explicit, it 
will be useful in the following to utilize the expression as 
an integral of multiple derivatives. Accordingly, we will 
define the integro-differential operator C'MSF(7r^(c;),p) 
which depends on the set E(Q), the ordering Tr^jg), and 
the limit p, and which acts on a set of variables ic in- 
dexed by the elements of E{Q). For notational sim- 
plicity, we again write it for E{Q) = {l,...,n} with 



'E(g) 



= vr' e Sn. 



OMSF(7r',p) 



)<---<^,'(„)<p 



i=l 



dL 



(H.16) 



We have 

E Omsf(^',p)/(^i 



,£„) =/(p,...,p), (11.17) 



for any nonsingular function / of n parameters, because 
the domain of integration becomes the cube £i < p, all i. 

Using eq. (jlLlSP in pi.lip . and the definition of 
CMSF(7r^(g),p), gives the desired expansion in terms of 
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topological graphs, 



OMSF(^£(g),p) E A 



(11.18) 



This is the final result of this section, which again should 
be compared with the percolation result (jll.4p . 

The expansion for MSF(p) that we have obtained is 
naturally organized as a low-density expansion, that is 
as an expansion in powers of p. For large graphs G, it 
becomes unwieldy, especially for p greater than around 
the percolation threshold (to the extent that a thresh- 
old can be associated with a finite graph, for example 
for a portion of a hypercubic lattice we can consider the 
threshold Pc of the infinite lattice). The corresponding 
expansion for the connectedness function Cx,y {p) in per- 
colation must produce the answer 1 when p —* 1, but in 
a very complicated way, as a sum of a large number of 
terms. For the MST, obtained from MSF(p) as p ^ 1, 
the probability that the path on the MST from x to y 
passes through z remains non-trivial in the limit, and is 
again given by a complicated set of terms. A general 
analysis of this sum on a large graph for p > pc would re- 
quire a resummation of terms to allow for the presence of 
the "giant cluster" in the corresponding percolation. In 
the Potts model formulation of percolation, this is done 
by giving an expectation value to the Potts spin. A for- 
mulation of such a resummation suitable for the MSF(p) 
problem will not be given in the present paper, which 
is consequently restricted to the region p < Pc on large 
lattices from here on. 



III. RENORMALIZED PERTURBATION 
EXPANSION FOR MSF PATHS 

In this section we describe how the exact low-density 
expansion obtained in the previous section may, for 
P Pc, be turned into a continuum theory (with a cut- 
off) to which renormalization-group methods may be ap- 
plied. This continuum theory is obtained from a naive- 
looking procedure of replacing the lattice paths repre- 
senting edges of a topological graph by continuum ran- 
dom walks, and neglecting the excluded volume require- 
ment that the graph be embedded in the lattice without 
using a vertex or edge more than once. The resulting per- 
turbation expansion is expected to be asymptotic rather 
than convergent. The expansion is closely related to that 
for percolation, so that the resulting Feynman diagrams 
and corresponding integrals can be compared with those 
of the latter. However, the expansion does not arise from 
the path integral of an action functional, so we do not 
technically have a "field theory of MSF paths," although 
we will show that many of the standard techniques of 
field theory remain applicable. In particular, we show in 



Section IIII El and Appendix [D] that the perturbation ex- 
pansion may be renormalized and RG methods applied. 
This enables us to calculate in Section IIIIFI the fractal 
dimension of paths on MSF(pc) as an asymptotic expan- 
sion in e = 6 — d, which we perform to leading order. 



A. The excluded volume constraint 

A lattice expansion in terms of topological graphs, 
such as (jll.4[) . is very close to describing a continuum 
theory. The only remaining roadblock lies in the sum 
over embeddings A, which carries an effective excluded 
volume constraint: edges of Q must be mapped to self- 
avoiding chains of lattice edges, and which must all be 
edge-disjoint (and hence also vertex-disjoint): no edge on 
the lattice may be used more than once. This is techni- 
cally more difficult to incorporate. If we drop this con- 
straint, we have a sum over "free embeddings" A, which 
map edges of Q to random walks in the lattice, which are 
allowed to intersect, and then Fourier transform tech- 
niques may be freely used. This is the starting point 
for a continuum theory: as is well known, the generat- 
ing function for random walks may be thought of as the 
propagator of a free scalar field. 

For percolation, the excluded volume constraint may 
be avoided by further modifications to the expansions 
given above. This is Essam's "p expansion" , given in 
Refs. [1^, [131 . There certainly seems to be no obstacle 
to extending the p expansion to the expansion for MSF 
paths derived in the following subsection, but we do not 
pursue this line of inquiry (which may be relevant for 
a mathematically rigorous reformulation of the results 
given here). Instead, in what follows, we assume we may 
drop the excluded volume constraint without difficulty 
or modification of our lattice expansion. This is because 
our ultimate aim is an RG calculation around dc — 
In the diagrammatic perturbation expansion, the impor- 
tant graphs for calculation of the exponents have ver- 
tices of degree 3 only. If we attempted to incorporate the 
excluded-volume constraints, perhaps following Essam's 
technique, additional diagrams (topological graphs) with 
vertices of degree four or more would enter, but these will 
be irrelevant close to six dimensions, as could be demon- 
strated by extending the techniques given below. 



B. Continuum perturbation expansion for MSF 
paths 

We first consider the percolation connectedness func- 
tions, defined via (jll.4p (and its generalization to n-point 
functions), which must be reproduced by the MSF(p) 
process. Because we neglect the excluded volume con- 
straint present in the sum over embeddings, for a given 
topological graph (which may now be referred to as a 
Feynman diagram) CJ, we may take the chains of edges 
produced by the lattice embedding to be independent 
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random walks. In the sum over such embeddings, it 
is natural to consider the Fourier transform with re- 
spect to the positions {i — 1, n), and to use 
the Fourier representation for the probability of a walk 
between two of the vertices; the latter takes the form 
l/(q^ + to) = Go{q_,to), the same as the propagator of 
a scalar field [23|. Here > is the "mass-squared" 
parameter, which depends monotonically on p; naively, 
to decreases to zero as p increases to Pc (however, this 
statement will be modified by perturbative corrections). 
Thus we make the substitution 

J2 PiiG<p]^iig,to), (III.l) 

where 

X Y[ {27:fS'' I (ke.t). - J2 -^--^-^ ) • 

Here J\f is the incidence matrix of Q under an arbitrary 
orientation of the edges: 

Afe^v = 1 if w is the head of e, 

= — 1 if w is the tail of e, and (HI. 3) 

= otherwise. 

The external momenta {koxt,i} are the Fourier conjugates 
of the positions {x^} of the graph's root vertices, and 
(kext)ti is the net external momentum flowing into vertex 
V. The momentum (i.e. wavevector) integrals are subject 
to a cutoff: each variable must obey \q_e\ < A. This 
cut off replaces the restriction of the integrals to a single 
Brillouin zone that is due to working on the lattice (note 
that in the latter case the propagators would be invariant 
under addition of a reciprocal lattice vector to any q^). 
Thus A is initially taken to be of order 1/a, where a is the 

spacing of the lattice points. Finally, a factor ^q^*'^'''"" 
has been inserted, to absorb other numerical factors that 
are omitted, and because this parameter will be renor- 
malized later. At this stage, go is strictly speaking of 
order one, but will be viewed as small in the perturba- 
tion expansion. At the same time, we will restrict the 
sum to topological graphs with vertices of degree three, 
except for the marked points x^, which are of degree one 
(thus we have go for each cubic vertex). Both of these 
simplifying assumptions can be justified because other 
contributions can be shown to be irrelevant (in a RG 
sense) near six dimensions, using the RG technology to 
be discussed later. 

Thus, the continuum expansion for the percolation 
connectedness functions becomes 

^{x.}(p)= E d{g)i{g,to)/A{g), (ih.4) 



where the sum is over diagrams with trivalent vertices 
and n "external" marked points on which d depends, 
and A{g) is the appropriate symmetry factor. For the 
function C^y{p), there is a similar expansion, in which 
we sum over graphs with two marked degree one exter- 
nal points at x, y, and with a single degree two ver- 
tex marked z (to replace e, as mentioned before); the 
symmetry factor becomes A'{g)^ and the power of go is 
T/((7)| — 3. In the same way that this function on the 
lattice was obtained by differentiating with respect to p, 
this function in the continuum can be obtained by dif- 
ferentiating the Feynman diagram expression for Cx,y(p) 
with respect to to (and no subsequent integration in the 
present case). In particular, this produces the correct 
symmetry factors. This operation gives the additional 
vertex with zero wavevector, but can be generalized to 
allow some momentum to enter at z, as given above. In 
field theory it is referred to as insertion of a mass- or 
(f>^ /2 operator, where ip would be the field corresponding 
to the degree one external points |29i| . 

Because of the equivalence of the d weights with those 
of the Q = 1 Potts model (at least for n — 1, 2), this 
perturbation expansion reproduces the standard one for 
the field theory of the Q — > 1 Potts model, which is 
usually obtained via the Hubbard-Stratonovich technique 
[23 , [23j . We emphasize that for the purposes of what 
follows we are prohibited from making any reference to 
an action functional due to the fact that the MSF path 
vertex cannot be expressed in terms of any local operator: 
instead we must phrase our argument entirely in terms 
of diagrammatic expansions. 

Now we turn to the path connectedness function 
Cx,y(p) for MSF(p), which is the probability that at pa- 
rameter value p there is a path on MSF from x to y 
passing through z. We treat the continuum version of the 
expansion in exactly the same way as we did for percola- 
tion, with the function y being the closest analogue. 
Then compared with percolation, we make the following 
modifications of the expansion: the diagrams are consid- 
ered for each ordering tt' of the costs of the edges e of 
the diagram, the d-weights are modified as they depend 
on the ordering through the requirement that the path on 
the MSF passes through z (and depend on marked points 
X, y, z), and the probability is modified to give the prob- 
ability for the ordering. Further, for the MSF, when we 
pass to the continuum perturbation expansion, the prob- 
ability for a given ordering of costs on the edges of the 
diagram must depend on squared-masses in place oi £e- 
These parameters are acted on by the integro-differential 
operator C'msf(7''',p), which finally sets all to p. Af- 
ter the change of variable to to, the operator C'msf(7''',p) 
may trivially be rewritten in terms of the {t^} variables: 
because every derivative is paired with an integral, the 
Jacobians involved in changing from the {i^} to the {t^} 
cancel. The only difference is that the ordering tt' ap- 
plies to the tf: in reverse: small corresponds to large 



9 



ig. Therefore, (|II.16|) becomes 

(III.5) 

The perturbation expansion for the MSF path vertex 
functions in the continuum is now 

E E {dMsAQW)/A'{G)) 

X C'MSF(7r',<o)/(a,{iJ). (III.6) 

The set of graphs involved are the same as those used 
in the differentiated two-point connectedness function of 
bond percolation (as described above), so we may com- 
pute this by starting with the expansion for y(to) and 
making the substitution, 

d{g)i{g,to) ^ 

E dMSF{QW)OMSYWM)I{Q,{te]) (ni.7) 

on a diagram-by-diagram basis. 

To summarize, the Feynman diagram rules for the 
MSF path correlation function, as specified in (|III.7p , are 
as follows: 

1) For each diagram Q G 5x.y,z contributing to the two- 
point correlation function between x and y of a cubic 
scalar field theory with a mass-insertion at z, we associate 
a mass-squared to < te < oo to each edge e G E{Q). 
For each ordering tt of these mass parameters, we act on 
the integrand with the operator Omsf{t^' ,to) defined in 

diim . 

2) After integrating over wavevectors, the contribution 
from each ordering is multiplied by the diagrammatic 
weight dMSF{QW)/ A'{Q), with c^msf defined in pi.Sp . 
and the sum of these over all mass parameter orderings 
is the contribution to the MSF path connectedness func- 
tion. 

Note that we must act with Omsf before any mo- 
mentum integrations are performed, since the latter may 
produce expressions that diverge at the upper limit of 
the integrations in Omsf- This situation could be reme- 
died by cutting off the domain of integration in Omsf 
to {^c} < A^, at the expense of complicating our RG 
calculation. 



C. Extraction of fractal dimensions and lowest 
order results 

The perturbation expansion that we have now obtained 
can be organized as a loop expansion: the lowest order 
contribution to y is order ©(gg); and is simply the 
diagram that takes the form of a path from x to z to 



y, which possesses no loops (cycles), while higher orders 
in go contain additional loops, one for each factor of g^. 
The lowest order result, in position space and at in = 0, 
takes the form (in this and the following, all separations 
like |x — y| are assumed large, 3> a) 

By contrast, the 2-point connectedness function at criti- 
cality, obtained as a single scalar propagator, is propor- 
tional in this order to l/|x — yj'^^^. [At zero-loop order, 
these results are the same for MSF(pc) and for critical 
percolation.] Dividing the two gives the conditional prob- 
ability that there is a path from x to y passing through 
z, given that there is a path from x to y. Integrating over 
z gives oc |x — yp. This is viewed as proportional to the 
total number of steps on the (lattice) path, even through 
the events of the path passing through the various z are 
not disjoint. The exponent 2 indicates that the fractal 
dimension of the path is 2, which is the correct result for 
a random walk. Thus we have shown that 

Dp = 2 (III.9) 

at zero-loop order. This will be found to be correct for 
d > 6, and also for d = 6 up to logarithmic correc- 
tions. The same dimension is believed to hold for paths 
on critical percolation clusters for d > 6, by similar field- 
theoretic arguments. Geometrically, it is because on large 
scales these clusters are trees, with no loops [^] , and hence 
are the same in the MSF(pc) process. 

In general, and specifically for d < 6 as we will show, 
the scaling exponents will be different. At Pc, quite gen- 
erally Cx y will have the scaling behavior 

C^bx,6y &-('^-2+,,)-(d-Dp)g,z^^ (III.IO) 

for any x, y, and z, and scalar 6, while the 2-point con- 
nectedness behaves as 

Thus these two functions determine two exponents rj and 
D-p for MSF(pc), and rj will be the same as for percolation, 
as we will explain shortly (77 = for d > 6). Then in the 
same way as at zero loops, we infer the fractal dimension 
D-p for the MSF path. From the geometric point of view, 
d — Dp is the co-dimension of the path. 

The exponent 77 describes the decay with distance r 
of the probability that two points are connected by a 
critical percolation cluster, namely ~ r~'-'^~^^^\ In a 
field-theoretic point of view, (d — 2 + ?7)/2 = a;^ is the 
dimension of the Potts field operator 0, while d—Dp = Xp 
is the scaling dimension for the path- vertex "operator". 

is related to the fractal dimension Dpcrc of the critical 
percolation clusters as the codimension x^ — d — Dpcrc, 

so i:»perc = (d + 2 - 77)/2. 

For the percolation function C^y, the corresponding 
0^ operator at e has dimension x^2 = d — Dgc, and Dgc is 
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the fractal dimension of the set of singly-connected edges 
on the path from x to y on the critical percolation cluster 
For d < 6 this set does not usually form a connected 
path. We see that the MSF path must include the singly 
connected edges, which leads to the inequality 

Dp > D,,. (III. 12) 

Because the function y is connected via differentiation 
with the change in connectivity with p (or to), the scal- 
ing dimension x^2 of the (jx^ insertion controls the length 
scale produced by taking p < pc, this length is the cor- 
relation length ^, and we can define the exponent v by 
C ~ (Pc — p)~'^ as p — > pc- It follows that Dgc = 
This discussion shows how the fractal dimension Dp , and 
others, can be extracted from the renormalized pertur- 
bation calculations. 

D. Beyond lowest order: breakdown of 
perturbation theory for d < 6 

The perturbation expansion for the MSF path connect- 
edness functions can be treated in a similar manner as 
that for standard field theories. A first step is to intro- 
duce one-particle irreducible (IPI) functions. A diagram 
is defined to be IPI if it does not become disconnected 
when a single edge is removed. Now for the MSF path 
connectedness function C^y, the (dominant) diagrams 
that contribute have a single edge emerging from x and 
y. For terms of order 0(50) (as 50 — > 0), the diagram 
possesses at least one loop (cycle). It can then be decom- 
posed into a chain of one or more disjoint IPI 2-point 
graphs, connected by single edges. The vertex labeled z 
is either inside one of the subdiagrams (subgraphs), or 
on one of the single edges. The IPI subdiagrams not 
containing z will be called self-energy diagrams. 

For the c?MSF-weight of such a diagram, it is easy to 
see that the weights associated with each IPI subdia- 
gram factor. This is because the MSF path must pass 
through each of the IPI subdiagrams in turn. For the 
self-energy diagrams, all paths through the subdiagram 
(which must be considered when evaluating the dMSF- 
weight) contribute a non-zero amount (all diagrams we 
consider are connected to x, y, z). Consequently, the 
factor in the dMSF-weight for the subdiagram is indepen- 
dent of the ordering tt'^ restricted to the subdiagram, 
and then the weight reduces to the same expression as 
in percolation. (This is not true, however, for the IPI 
subdiagram that contains the path vertex at z.) The ap- 
plication of the Omsf operator and the sum over order- 
ings tt'^ can now be carried out using pi.l7[) . Then the 
contribution of such a self-energy diagram is the same as 
in percolation. The self-energy diagrams can be formally 
summed to all orders in perturbation theory to yield the 
self-energy S(q, to), and then each of the two series of al- 
ternating Go(q, to)'s and S(q, to)'s can be summed as a 
geometric series, giving the full Green's function G(q, to), 

G(q, to)-' = Go(q, to)-' - S(q, to) (111.13) 



(Dyson's equation). We pause to point out that the 2- 
point connectedness function for MSF(p), in which we do 
not require the path on the MSF to pass through any par- 
ticular point z, is similarly shown in this diagrammatic 
point of view to be the same as in percolation, and is 
given by G(q, to). Consequently, the exponent rj defined 
above for p = pc must be the same as the similarly- 
defined exponent in percolation. 

The MSF path connectedness function, with zero mo- 
mentum entering at the path vertex, can now be written 
formally as 

jd''zCly{t0) = 

J -^^'^^^^^''^ ^^"1' ^0)^^'^'^''^ (q^ 0; ^o)G(q, to), 

(III. 14) 

where F^^'^^^ (q, 0; to), which we call the path vertex 
function, is the Fourier transform of the sum of all IPI 
diagrams with two external points (connected to x, y), 
plus the path vertex at z, which has here been assigned 
zero momentum. (The generalization to r'^'^^(qi, q2; to) 
should be obvious.) Diagrams contributing to F*^^'^^' 
are depicted in figure [TJ Similarly, we also define, for 
the iV-point connectedness functions G'^^(qi, . . . , qA?) 
without the path vertex (TV = 2, 3; G^^^ = G), and the 
2-point connectedness function with a mass (or ^^) inser- 
tion, which are the Fourier transforms of the percolation 
functions G^^) = Gx,y, G^^-i) = G^^y: 

TV 

GW({q.};io) = n^(q-^o)-r^'^^(H};^o); 

i=l 

G(^'i)({qa;q;to) = 

N 

nG(q„to)-r(^'l)({qa;q;to). (III.16) 

i=l 

In these functions, a 5-function that sets the total 
wavevector to zero has been removed, and {qi} stands 
for the ordered set qi , . . . , qAr . This causes a minor dif- 
ference in notation from that for the path vertex function 
above: in the functions G or F containing or + 1 
wavevector arguments, one of the wavevectors could be 
eliminated, which is what was done in r^^'^^^^ above, and 
we occasionally do this for the others also without fur- 
ther comment. The functions F^^' and F^^'^^ are called 
the IPI vertex functions (of the types indicated). We 
identify F(2)(qi, -q^; io) = G(qi;to)-i. 

The problem of calculating the path exponent Dp has 
now been reduced to the calculation of the IPI path ver- 
tex function F*^^'^^^ . The external lines G are the same as 
in percolation, because of the factorization and ordering 
independence of the weights for the self-energy diagrams. 
(Similar, but more subtle, factorizations play an impor- 
tant role in the later part of the argument also.) The 
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FIG. 1: Perturbation expansion for the path vertex function 
given in (|III.6|I . We represent the path vertex by the open 
circle connected to a wavy hne. Note that these diagrams do 
not include any notation corresponding to the operation of 
Omsf- The lowest non-trivial term is the second one on the 
right-hand side. 



path vertex, on the other hand, is not the same as the 
mass-insertion vertex in percolation which it resembles. 

To illustrate the perturbation expansion, let us now 
evaluate the first correction, of order gg, to F^^'^^^ (see 
Fig.[T]) at zero external momentum. From the rules given 
above, this correction comes from the graph with three 
propagators connected to form a triangle, and the con- 
tribution is (note that A' = 1 for this graph) 

^MSF(A,to) = 

dMSF(A|^)OMSF(^,io)/A(A,{ii,t2,t3}), (ULIT) 

where the last integral is 



/A(A,{ii,t2,t3}) = 



(2^)^ (k2-f ti)(k2+i2)(k2+<3)' 



(in.18) 



Here is the mass-squared on the edge not adjacent to 
the path vertex. Our momentum-space rules required us 
to integrate up to radius A. Since the diagram is eval- 
uated at zero external momenta, it is symmetric under 
permutations of the {t^} and 

OMSF{n,to)lAiA,{h,t2,t3}) = ^/A(A,{to,to,to}). 

(in.19) 

Applying (jn.8[) . we have (iMSF(A|7r) = for the two or- 
derings in which ti and t2 > is, and (iMSF(A|7r) — —1 for 
the other four orderings. Thus the result is 



/MSF(A,to) = --^lA{A,{to,to,to}). 



(in.20) 



For the corresponding simple mass-insertion ver- 
tex in percolation, the result would be instead 

-lAiA,{to,to,to}). 

Dropping and numerical factors, the contribution to 
the path vertex function behaves like the integral 



I'a 



A i,d-l 



{k^+W 



(in.21) 



At present, we are interested in this for fixed A as to 0, 
so as to reach p = Pc- We see that for d > 6, (|in.2ip be- 
haves as A^~^ as A — > oo, and is finite (for any A) as 



to ^ 0. But for d < 6, the reverse is the case: the in- 
tegral converges as A — > oo, but diverges (for any A) as 

tff ^''^^ as to 0. In the borderline case d = 6 the inte- 
gral diverges logarithmically at both ends. For non-zero 
external wavevectors, the dependence of the integral on 
A is the same in all cases. Note that similar statements 
apply for percolation; only the numerical prefactor is dif- 
ferent. 

There are similar results for diagrammatic contribu- 
tions with more loops. Simply counting the number of 
propagators and integrations gives the "superficial de- 
gree of divergence" , which for the path vertex function 
is always ~ gpA'*"^ or Sq^q'^"^''^^ for d > 6 and d < 6 
respectively, raised to the power of the number of loops 
(independent cycles) in the diagram, as in the one-loop 
example above. Note that this is the same as if the Omsf 
operator were absent, because Omsf leaves the overall 
degree (in k, at to ~ 0) of the integrand unchanged. The 
consequence for perturbation theory at fixed A and d > 6 
is simple: each term in the perturbation expansion of F^^ 
for MSF(p) {p ~ Pc) is finite as to 0. This is true for 
the self-energy diagrams on the external lines also, and 
the value of I](q, to) at q = determines an effective shift 
in the value of to that corresponds to p = Pc' Pc must cor- 
respond to the value of to such that to — S(0, to) = 0, and 
there are also other finite changes in the normalization 
of G. (In this case one would wish to sum up self-energy 
insertions in the lines inside of F^^ also. However we 
have not shown that these take the same form as on the 
external lines. This will be addressed below.) But the 
consequence is that in each order, F^^(q, to) ^ 0(1) as 
q — > at p = Pc, while G'(q, to)~^ oc q^. This in turn im- 
plies that there is no change in the exponents from their 
lowest order values, rj = and Dp = Dgc = 2. Note that 
here we disregard the possibility that the sum of an in- 
finite number of finite terms might diverge, which might 
invalidate the conclusion. 

For d < 6, this perturbative argument breaks down 
as the corrections become arbitrarily large as to ~^ 0; 
in particular in the region to < .^o^'^ '^^ (the Ginzburg 
criterion). In order to handle this, the use of RG tech- 
niques becomes essential. These techniques effectively 
re-sum and redefine the expansion. There are several 
formulations of the RG. These may be divided into two 
classes. One class of particularly powerful techniques is 
the field-theorists' RG, in which the aim initially is to 
take A — > oo (or a — > 0) at fixed separations or mo- 
menta, in such a way that the limits of the correlation 
(or connectedness) functions exist, thus recovering a true 
continuum theory. This is called renormalization of the 
theory. Subsequently, the renormalized theory is used 
to set up the RG, and calculate exponents for d < dc- 
The leading alternative is the Wilsonian RG, in which 
the cutoff is kept finite. The Wilsonian RG is more diffi- 
cult to use for higher numbers of loops. Both approaches 
lead to equivalent results for physical quantities such as 
exponents for d < dc, where dc = 6 for percolation and 
MSF(pc)- In this paper we will follow the approach of 
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the field theorists. 

For d > 6, we can see from above that the effective 
expansion parameter is (^gA'*"^. As A ^ cxd, it is then 
necessary to make go ^ such that ggA''"^ does not 
diverge. In fact the situation is even worse than this 
would suggest: there simply is no rational way to define 
the limit A ^ oo so that the connectedness functions at 
fixed q and (for example) p = Pc have finite limits, with- 
out introducing an infinite number of parameters. This 
is referred to as non-renormalizability of the perturbation 
expansion. But by keeping the cutoff A finite, and using 
the Wilsonian point of view, we can see that the expo- 
nents in this region take their simple zero-loop values, as 
indicated above. Accordingly, we concentrate on d < 6 
from here on in this article. 



E. Renormalizability of the theory 

In this section we outline our proof that the MSF path 
vertex theory for d < 6 may be consistently renormalizcd. 
The full technical details are in Appendix iDl In the in- 
terests of making this article more accessible to readers 
without a field theory background, we take a somewhat 
pedagogical approach in discussing the renormalizability 
and the RG calculation in the remainder of this paper. 
Of course, we do not have the space here for a full de- 
scription; the interested reader is directed towards any of 
the standard textbooks such as [13, HH . 

We saw above that the Feynman integrals associated 
to diagrams for certain vertex functions are superficially 
divergent as A ^ cx) in six dimensions. In fact, closer in- 
spection reveals that subintegrals (integrals over a subset 
of the loop momenta k, holding the others fixed) may also 
be superficially divergent, and this can occur even when 
in integrals that are superficially convergent as a whole, 
showing that they do not converge after all. However, 
the superficially divergent integrals (or subintegrals) are 
associated only with (sub-)diagrams that are, topolog- 
ically at any rate, of the form of the vertex functions 
r(2), r(3), r(2.PV)^ or r^^.i). These correspond respec- 
tively to the self-energy, cubic coupling, MSF path, and 
mass-insertion vertex functions. The first two of these 
suggest a possible way to eliminate the divergences: add 
the divergent terms to the mass-squared to and the cou- 
pling go, respectively, and define renormalizcd quantities 
t and and then insist that these are the "physical" or 
measurable parameters at long length scales. There is 
also a subleading divergence in the self-energy of order 

, which perhaps can be removed similarly by rescaling 
the field, and hence the Green function (this effect also 
enters the definition of t). The divergences in F^^'^^^ and 
r^^'^^ can be handled similarly. 

This procedure works for conventional field theories at 
their critical dimension; one such case is the theory of 
percolation at six dimensions. It is important to recog- 
nize why it can work. That is because every occurrence 
of, for example, a self-energy subdiagram within another 



diagram occurs with a d weight that can be factored as 
the d weight for the subdiagram, times that for the quo- 
tient diagram, in which the subdiagram is contracted to a 
single vertex. Further the integrations over wavevectors 
in a Feynman integral have the property that the integral 
for a subdiagram always has the same form, independent 
of the larger diagram of which it is a part. Then the 
subintegral for any self-energy subdiagram has precisely 
the same divergence wherever it occurs, independent of 
the larger diagram of which it is a part. Meanwhile 
the quotient diagram has the form of a lower order dia- 
gram. This enables us to write the leading divergence as 
a correction to the bare mass-squared to that is context- 
independent, and therefore meaningful. The same has to 
be true for the other divergent subdiagrams (or "renor- 
malization parts"). For percolation, the factorization of 
the d weights can be easily seen in the Q ^ 1 Potts model 
formulation, in which the d-weights arise from contract- 
ing together tensors, and then the factorization for sub- 
diagrams is automatic (and similarly for other local field 
theories). It is not immediately obvious that this will 
hold for our MSF theory, because: (i) The d weights are 
replaced by dusF weights, which depend on an ordering 
TT, may not factor in the fashion required, and in fact 
for some orderings do not factor; (ii) The Omsf operator 
and summation over orderings raise similar questions. 

In Appendix [D] we undertake a careful study of these 
questions. We find that the degree of divergence of a 
subdiagram for a renormalization part is the same as it 
would be for the corresponding subdiagram in percola- 
tion for certain orderings, and for these the dMSF-weight 
exhibits the desired factorization properties. Indeed, for 
the self-energy and cubic coupling renormalization parts, 
the divergent part has exactly the same coefficient as for 
percolation. This holds also for part of the subleading 
divergence in the self-energy case, but there is also an- 
other subleading part in that case which does not have 
these properties. That part is problematic, as the pro- 
gram above provides no apparent way to remove these 
divergences. However, we eventually find that all such 
terms cancel, not for a single diagram, but in the sum 
of diagrams of a given order. In the remainder of this 
discussion, we will take that for granted, and so continue 
as if there are no such divergences. 

These observations then allow us to absorb all the di- 
vergences into the quantities mentioned above. More for- 
mally, this amounts to subtracting off the superficially- 
divergent contribution for each subdiagram of a diagram 
(including that, if any, for the diagram as a whole). Af- 
ter doing so, we should prove that the remaining integrals 
are actually finite. Here again, we cannot simply appeal 
to the usual field theories, as we have modified the Feyn- 
man integrands, and so the proofs must be reconsidered. 
We complete the proof usin g th e Schwinger parametric 
integral formulation 



32 , |33|, l3J, l35|, l36(] , and a theorem 
37[. This then completes the proof 



by Bergere and Lam 
of renormalizability of our perturbation expansion. 
The renormalization procedure removes the divergent 
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parts of the original Feynman integrals. It does not 
uniquely fix a finite part that may also be subtracted. 
This part may be determined by giving some renor- 
malization conditions obeyed by the renormalized vertex 
functions. A convenient choice for the following is to de- 
fine the values (and a first derivative) of these functions 
at zero renormalized mass-squared, t = (corresponding 
to p = Pc), and a non-zero wavevector of magnitude k. 
For dimensional reasons, one or other of t and k must en- 
ter. However, we also mention the scheme of dimensional 
regularization and minimal subtraction, in which such 
renormalization conditions are not used. We adopt the 
present scheme in order to keep things relatively trans- 
parent. 

The RG is now introduced by obtaining an equation, 
the RG equation, describing how the renormalized vertex 
functions behave under a change in k. As k decreases, the 
effective coupling g changes, and may reach a non-zero 
K-independent fixed point. This is then used to calculate 
the exponents for scale-covariant behavior of the vertex 
functions or correlation functions. The fixed point is at 
of order e — 6 — d, and so the expansion in powers 
of g is traded for one in powers of e. This expansion is 
essential to obtain useful finite results for d < dc = 6- In 
this way we will obtain the exponents to order e via a 
one-loop RG calculation. 



of corresponding wavevectors.) We will also now intro- 
duce dimensionless versions of the bare and renormalized 
couphngs go, g: 



(47r)'^/2 



g^K ^ 
(47r)'^/2 



(III.24) 



We introduced an angular factor (47r)'^/^ in the above def- 
initions for later convenience to simplify expressions. For 
the vertex functions not containing the path vertex, the 
functions and their renormalization is exactly as perco- 
lation, and this is also true of the following calculations; 
we include some details anyway to provide checks on the 
calculation. 

The RG equations are obtained from the observation 
that the bare functions Fq are independent of k when 
written in terms of go, so ndVo/dn = at fixed go, to, A. 
Using the definition of we obtain 



d_ 

dn 



I3{u] 



d N 



= 0, 



(III.25) 



and 



F. RG analysis at one-loop order 

In the preceding sections and appendix [Dl we have 
proved that the diagrammatic expansion for the MSF 
path theory is renormalizable. These proofs were techni- 
cal, but having established this fact, we are free to make 
use of standard RG methods such as those discussed in 
Ref. [3l|. We continue to take a rather pedagogical ap- 
proach in this section. 

As we explained in sections IIII Dl and |D1 renormaliz- 
ability of a theory means that we may absorb the strong 
A-dependence of all correlation functions into a finite 
number of parameters and the overall scale of the cor- 
relations, at the cost of introducing another scale k. In 
the scheme we use, in which the renormalization condi- 
tions are at zero renormalized mass-squared t = and 
non-zero wavevector of order k, the precise statement is 
that functions g, Z^, Z^2, Zpy exist such that (here we 
append subscripts to denote the "bare" vertex functions 
as constructed above, with cutoff A) 



Z;'/'(5o, A)Z^. (50, «, A)F^"^"^({q,}, {q,}; go, to, A) 



(NX), 



= F 



({qj,{qi};5,«:). 



(III.22) 



(2, PV) 

^^(ffo, «, A)Zpv(go, K, A)F^ ' '{{q_i};go,to,A) 

= T'^^:''^\{q,};g,n), (111.23) 

where Tn are independent of A as A oo with g, t — 
fixed, up to corrections vanishing in this limit. (Here L 
is the number of insertions of cf)^, and {q_j} is the set 



xF(,^'P^)({qJ,H,«:) 
= 0, 



(IH.26) 



In each of these equations the first and second partial 
derivatives are at fixed u and fixed n, respectively. The 
RG (3 and 7 functions appearing in equations (|III.25|) 
(IIII.26P are defined as 



70 (w) 
(u) 
7pv(m) 



du 
d log Z^ 



so. A 

dlogZx'i 



dn 
d log ZpY 



Ok 



So.A 



So, A 



(IH.27) 
(in.28) 
(IH.29) 
(IH.30) 



and are finite as A ^ cx) 31]. Hence in the limit they are 
independent of A, and so also of k, because u, (3, and all 
7's are dimensionless; they are simply power series in u. 

We then impose the following renormalization condi- 
tions, which are those we reached in Appendix [Dl but 
written now with t = Q and u in place of g. These serve 
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to fix the dependence of the parameters on one another: 



d_ 



R 

.(2) 



|q|=K 



t = 0, 

1, 





k) 


SP 


= 3. 








n) 


SP 


= 1, 








n) 


SP 


= 1. 



(III.31) 




FIG. 2: The IPI one-loop diagrams for vertex functions that 
are ultraviolet divergent in six dimensions: a) self-energy, 
b) cubic interaction vertex function, c) mass-insertion vertex 
function, d) MSF path vertex function. 



Here SP denotes a symmetry point of the external mo- 
menta qi, ^2^ qs, at which |qi| defined in Ap- 
pendix |D] (though the precise definition is unimportant). 
These conditions are now used to determine /3 and the 
7's from the perturbation theory expansion in go of the 
IPI vertex functions Tq with a fixed cutoff A. The ex- 
pressions make sense provided 50 is sufficiently small. We 
will calculate to one-loop order, which means that only 
the one-loop diagrams for the renormalization parts need 
to be calculated. This will give results for exponents to 
first order ine = 6 — (i>0. (More generally, computing 
to 0{e'~) in the e-expansion requires computing all the 
renormalization parts with C or fewer loops.) 

Then the instances of equation pil.23p with which we 
need to deal are, to 0(50): 

rg^ (q, u, k) = ((q^ + to) - S(q)) , 
S(q,.go) = d2glI2{^)\ 



SP 



^3/2 / 

[90 



-r(24)/ N 



pPV 



{u,k)\ 



SP 




d29lh 



SP/ 



^ c^MSF(Al7r)C'MSF(7r)/3(7r)|sp j . 

(III.32) 



Here ^2 , ^3 are the percolation d- weights for these one- 
loop diagrams, d2 = —1, d^ = —2, and dpv will be evalu- 
ated in a moment from the weights d{Q) and the operator 
^msfCmsf; the values will be substituted only at the end 
of the calculation. The negative sign in the equation for 
Fp^ arises because of Dyson's equation. 



In equation (jIII.32p . I2 and are Feynman loop in- 
tegrals which we now evaluate. To leading non-trivial 
order, we can evaluate them at six dimensions, retaining 
only the terms that diverge quadratically or logarithmi- 
cally as A — > 00. These terms may be extracted by any 
of the standard techniques for evaluating Feynman inte- 
grals, including those in Appendix iDl we refer the reader 
to [3il.l38l.l39l in particular. Denoting this approximation 



by ~, we find 



^2(q) 



d'^'k 



1 
2 
1 

2 (47r)''/2 



1 



{2ttY k2(k-Kq 

1 / .2 



A^ 



q- A 
Trlog — 
6 q- 



(IH.33) 



The factor of 1/2 appearing in /2(q) is a diagrammatic 
symmetry factor, = 2 (in the other diagrams, A or 
A' — 1). Note that the bare propagators appearing in 
these integrals should have mass-squared to = ic deter- 
mined so that t = 0, but because tc is 0(50) in^-Y 
consistently neglect its presence here. The other integral 
is 



'3lsp 



d'^'k 



1 



{2ttY k2(k + qi)2(k-q2)2 
1 1 



SP 



2 {AttY/- 



A2 

■log—. 



(HI. 34) 



For the path vertex function, as we saw above, 
C'MSF(7r, tc)^3|sp is independent of the ordering tt. (Un- 
fortunately, this property does not hold to higher orders 
in perturbation theory for F^^.) The result for a single 
ordering is 



Cmsf(7I', tc)^3|sp = ■^^3|sp- 
Then for brevity we define dpy such that 



(IH.35) 



^ dMSF(A|7r)OMSF(7r, tc)-?"3lsp = c;pv^3lsp- (HI. 36) 

Hence dpv — \ X^tt ^MSF(A|7r) — —2/3, as we saw above. 

We may now solve by requiring that the renormalized 
proper vertices defined on the right-hand side of (1111.32^ 
satisfy the normalization conditions (IIII.31|) up to terms 
of 0{uq). We obtain to this order 



tc = |c^2^^oA^ 



Zpv 



= 1 - ^^2^0 lo; 



{d3 
A 



^d2)ul \oi 



1 + (5^2 - dpY)ul lof 



A 



(HI. 37) 

(IH.38) 

(IH.39) 
(IH.40) 
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Finally, to this order Z^i is the same as Zpv except that 
di replaces dpy- Note that the Z's are functions only of 
the dimensionless variables mq and k/A, and we can set 
d = 6 (so e = and — Ug/(47r)'^)) in the one- loop 
terms, but not in zero-loop terms. 
The RG functions are 

+ C'(u^e^i3), (III.41) 
70(u) = id2^i^ + 0(^i^eM2), (III.42) 



7Pv(u) 



Ida - dpv) + C(^^^ eu^). (III.43) 



The fixed points of the RG are the values of u = u* at 
which /3(m) = 0. Clearly, one fixed point is at u = 0, but 
is unstable to the introduction of the cubic coupling for 
d < 6, because u grows as k decreases, corresponding to 
the behavior at larger length scales. In six dimensions, u 
approaches zero logarithmically as k decreases, because 
the coefficient oi is da /4 — da — 7/4 which is positive. 
Below six dimensions, there is another zero of /3 which 
results from the competition between the two terms in (3, 
at 

d2 - 4d3 

Note that (u*)^ is positive for d < 6 (e > 0). For d > 6, 
this fixed point is not relevant to percolation or MSTs. 

The values of the 7's at the fixed point value of u give 
the "anomalous dimensions" of the various operators (ex- 
cept that in the case of 0, the anomalous dimension is 
70/2). These are the difference of the total dimensions x 
of the operators from their canonical or engineering di- 
mensions, which are the zero-loop values discussed earlier 
[3l| . The co-dimension d — x gives the fractal dimension 
of the associated geometric object (set of points). The 
most interesting dimension for us is that of the path ver- 
tex, which is xpy — d — 2 + 7pv. The codimension yields 
the fractal dimension Dp — 2 — 7pv of the path on the 
MSF(pc), as discussed in section HIl CI We find 

(d2 — 6dpv)e 



7pv(m*) 



that is 



3(d2 - 4d3) 



+ 0{e^), (III.45) 



(III.46) 



This is the main quantitative result of this paper. 

The other 7's produce fractal dimensions related to 
properties of percolation, which also apply to MSF(p). 
First, 



70(m*) = 77 



doe 



3(d2 - 4d3) 



0[e 



(III.47) 



Hence for the fractal dimension -Dpcrc of the critical per- 
colation clusters, which is D^r 



d — X0 , we find 



(d-f 2-?7)/2 



lOe 
'21 



(III.48) 



The other 7 is 7^2, for which the value at the fixed 
point can be obtained from the formula for 7pv by re- 
placing dpv by d2, that is 



702(u*) = 2-Ac = |^ + O(e' 



(III.49) 



and so for the fractal dimension of the set of singly- 
connected edges (see section ITlI Cp we find 



i^sc = jy-^ =2- 



5e 
21 



(III.50) 



The values we have obtained for both exponents rj and v 
a gree with those in the literature on percolation, to order 
£ |29l. [40I I4H ■ which provides a check on our calculation. 

The comparison of Dp with Dsc raises some questions 
of inequalities obeyed by Dp. There are also some other 
fractal dimensions defined for paths on critical percola- 
tion clusters which have been studied. These include 
Dinim the fractal dimension of the shortest path on the 
cluster between the given points, and Umax, the frac- 
tal dimension of the longest (self-avoiding) path between 
them [?! . Here the length of the path is the number of 
edges of the lattice that it traverses. Then the inequali- 
ties are fairly obvious: first, because all these paths must 
pass through the singly-connected edges, Dgc is the small- 
est of all, and the remaining inequalities 



(IH.51) 



follow from the definitions. To order e, one has -Dmin = 
2 — e/6 and Dmax = 2 — e/42 [i^l, and all the inequalities 
are obeyed strictly by the results to this order. Dp is 
close but not equal to -Dmin- 



IV. CONCLUSION 

The results of this paper fall into three main parts. 
First, we constructed an exact expansion for the Kruskal 
process, or spanning forest MSF(p), in a series in powers 
of p, which terminates for a finite graph, and is anal- 
ogous to a low-density expansion for percolation, or a 
high-temperature expansion in a statistical mechanical 
model. The expansion is for the probability that the 
path on the MSF from x to y passes through a vertex 
z. Second, this expansion was used to obtain a contin- 
uum formulation (with cut-off) for p < Pc (where Pc is 
the percolation threshold) in terms of Feynman diagrams 
(the region p > pc presents additional technical prob- 
lems, and we will not discuss these further here). This 
expansion was then shown to be renormalizable to all or- 
ders in perturbation theory, so that the limit of infinite 
momentum-space cutoff (or zero lattice spacing) can be 
taken. Third, the renormalized perturbation expansion 
was used to calculate the fractal dimension of any path 
on MSF(p) at p — Pc, to first order in e = 6 — d, for 
d < 6: Dp = 2-e/7 + 0{e'^). For d > 6, Lip = 2. If the 
"superhighways" idea is correct, then the same Dp also 
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applies to the region p > Pc, in which we expect the path 
dimension to be independent of p on large enough length 
scales. 

It is important to realize that it was by no means ob- 
vious at the outset that such a field-theoretic renormal- 
ization process would be possible. The problem is not 
obviously given by a local field theory, and our expan- 
sion is not based on an action principle (at least, not 
in any apparent way). Optimization is generally a non- 
local process as it involves making comparisons among 
(sums of) costs globally; however, this is also true when 
one wishes to minimize a Hamiltonian, even if its param- 
eters (corresponding to costs) multiply local interaction 
terms. For minimum spanning trees, the definition of the 
allowed or "feasible" configurations (i.e. spanning trees) 
is not local either. It was not obvious that the expansion 
would be renormalizable like that of a local field theory. 
Indeed, in the end our procedure worked thanks to un- 
expected and non-local cancellations of some subleading 
divergences (see Appendix lD 6p , for which we are unaware 
of any analogs in local field theories. Undoubtedly the 
underlying reasons for this success with MSTs should be 
found in the applicability of Kruskal's greedy algorithm 
and its connection with percolation. 

The calculations can be extended in various ways. The 
exponents can be calculated to higher orders in e, with 
increasing effort required for each additional order. The 
path vertex function, and not only its scaling dimension, 
can in principle also be studied, as can more general cor- 
relation functions with path vertices and mass-insertion 
vertices. In six dimensions, there are logarithmic correc- 
tions to the simple scaling with Dp ~ 2 that holds for 
dimensions bigger than six, and these arc calculable. 

Independently of these applications of the renornial- 
ized perturbation expansion, the exact lattice low-density 
expansion could be studied in low orders (say, the first 
thirty terms) in any dimension d, as is conventionally 
done with high-temperature series. This would provide 
another way to obtain scaling dimensions for correlation 
functions. Such techniques are frequently very accurate. 

A further question is the Borel summability of the per- 
turbation expansion, or of the e expansion for the expo- 
nents. If an asymptotic expansion of a function is Borel 
summable, then it uniquely determines that function. If 
a few terms of the expansion are available, and it is be- 
lieved to be Borel summable, then an improved estimate 
for the quantity of interest, such as an exponent, for a 
non-zero value of the parameter (say e = 1) can be made, 
and for critical exponents these values may be very accu- 
rate (comparable with high-temperature series methods) . 
For percolation at threshold, the asymptotic high-order 
behavior of the perturbation expansion has been shown 
in Ref. [43| to have the form that is a necessary con- 
dition for the expansion to be Borel summable. These 
results also apply to our theory, but again we also need 
a similar result for the path vertex function. It would be 
interesting to find a technique to estimate the high order 
behavior of our expansion. 



In conclusion, the introduction of the Kruskal process 
and geometric object MSF(p), based on an optimization 
problem, provides a rich area for study not unlike conven- 
tional critical phenomena. At p — pc, many techniques 
can be applied to it. It illuminates numerical work on 
such problems as optimal paths and transport in random 
media. 
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APPENDIX A: LOW-DENSITY EXPANSION IN 
PERCOLATION 

This Appendix summarizes various results concerning 
the low-density diagrammatic expansion for percolation. 
Some of these results have appeared elsewhere in the lit- 
erature (25l . [26j . but we reproduce them here in order 
to introduce notation and terminology, and because the 
derivation in Appendix[C]below closely follows that given 
here. 

In section lA II we define the low-density graphical ex- 
pansion for two-point connectendess functions, which is 
generalized to the case of n-point functions in section 
lA 21 These sections summarize results given in [2^ [2^ . 
Finally, in section [A 31 we prove that the expansion is the 
same as that obtained from the conventional description 
of bond percolation via the low density (high temper- 
ature) series for the Q-state Potts model, in the limit 
Q — > 1. The principal results are the graphical expan- 
sions (jA.9p . (|A.16p with diagrammatic weights given in 
simplest form in (jA.Up . (jA.lSp . 



1. Essam's construction 

Essam's expansion for percolation [2^ [2^ is based on 
the principle of inclusion and exclusion from combina- 
torics 44] . As this may not be familiar to all readers we 
summarize it here. We start with a set of events {Xi} 
indexed by i in an index set /. In order to calculate 
probabilities later we introduce the indicator function 



1 Xi true, 
Xi false. 



(A.l) 



The principle of inclusion-exclusion is the expansion 



iel -' iel i,jel:i<j 



^ E (-i)i^'i+ii A X 



(A.2) 
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where V denotes a logical OR and A denotes a logical 
AND. An analogous series may be obtained for the con- 
junction of all the events, by using De Morgan's law 
^(V.,Xj) = /\i{^Xi), where -. denotes logical NOT. This 
yields 



be rewritten 



I 



iei 



= 1 -: 



iei 



I'CI 



(A.3) 
(A.4) 



Note that in this case /' may be the empty set, so the 
first term of this series is 1. 

We apply this to bond percolation at a parameter value 
p by first investigating the two-point connectedness func- 
tion, defined as 



Cx,y(p) = (Ic(x,y)|p). 



(A.5) 



where Ic(x, y)|p stands for Ic(x,y)|p — 
I [x, y connected by edges of cost < p] and the an- 
gle brackets denote an average with respect to all 
realizations of the edge costs. Defining Fx.y to be the 
set of all self-avoiding walks on the lattice between x 
and y, we may write 



Ie(x,y)|p =1 
where we define the event 



V (^^p) 



il <p) = max 4 

V ee7 



< 



(A.6) 



(A.7) 



i.e, we require all edges e on the path 7 to be present by 
the time the parameter is raised to the value p. Using 
equation (|A.2p to expand the right-hand side of (|A.6[) by 
inclusion-exclusion yields 



Ic(x,y)|p= Yl 



(_l)|r'|+ii 



7'er' 



(A.; 



We obtain an expansion in terms of graphs from (jA.8|) 
by grouping together all terms that test the same set 
of edges on the lattice; the terms in the series are now 
indexed by graphs G, each obtained as the union of some 
set of paths T' C rx,y (possibly from more than one such 
r'). We say that such a set F' covers (the edges of) G. 
Because the paths are self-avoiding walks, all the graphs 
G generated from such unions must be vertex-irreducible: 
removing any vertex from the graph must leave at least 
one of the points x, y in each connected component. Let 
the set of all such graphs with the marked vertices x, y 
be Gx,y 

Letting rx,y(G) for G G Gx,y denote the set of paths 
on G connecting the root points x, y, equation (jA.Sp can 



=(x,y)|p = 

H (-l)l^'l+il[r' covers G]l[G < p]. 



GeG^.y r'cr^,y(G) 



(A.9) 



The average over the costs can be performed immedi- 
ately. Referring back to the definition IjA.Sp . we obtain 
the graphical expansion 



Cx,y(p)= E d{G)PT[G<p], 



(A.IO) 



GeGx 



by introducing 



d{G e Gx,y) = (-l)'^''+^I[r' covers G], 

r'cr^,y(G) 

(A.ll) 

which is independent of the parameter p, and is referred 
to as the d-weight of the graph G. This expression for the 
d- weight can be seen to possess the topological invariance 
property mentioned in section III Al the insertion of any 
number of vertices of degree two (or, equivalently, replac- 
ing edges of G with paths of edges) does not change the 
set rx,y(G) of paths connecting the root points, or the 
value of I[r' covers G] for any of the subsets F'. 

The definition of d{G) may be extended to cover the 
case where G is any two-rooted graph as follows: if G con- 
sists of more than one connected component, there is no 
way to cover all its edges with paths connecting the roots, 
so d{G) — 0. Note that, because the covering criterion is 
defined in terms of the edge set only, addition of isolated 
vertices does not change a graph's d-weight. Similarly, if 
G is not vertex-irreducible, by definition some edges — 
the "tadpoles" or "dangling ends" — cannot be covered 
by a self-avoiding path, since backtracking is forbidden, 
so again d{G) — 0. Since d{G) vanishes for these addi- 
tional cases, the sum in (jA.lOp may be extended to all 
two-rooted subgraphs of the underlying lattice. 

We may make further progress if we remark that the 
preceding derivation also applies to connectedness func- 
tions on an arbitrary graph A instead of the whole lattice; 
the sum in (jA.lip is then over appropriate subgraphs of 
A. We denote this connectedness function by Gx.y(A,p). 
Equation (jA.lOp generalizes to 

Gx,y(A,p)= J2 d{GE')Pr[GE' <p], (A.12) 

E'CE{A) 

where Ge' is the subgraph of A consisting of all vertices 
of A and a subset E' (- E = E{A) of its edges. Evaluating 
(1X12)1 at p = 1 yields 

I[E{A) connects x,y] = Gx,y(A, 1) = ^ d{GE'). 

E'CE{A) 

(A.13) 
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In the definition of d{G) ior G = G e' , defined by a subset 
E' of tlie edges of A, we may note tliat vertices of A 
incident on no edges can be deleted without changing 
d{G). Now because A is an arbitrary graph, and the sum 
in (|A.12[) is over over all subsets of E(A), we may easily 
invert this sum by Mobius inversion [4j|, which for the 
present case is related to inclusion-exclusion. We obtain 



d(G) = 



^ h[E' connects x, y] 



E'CE(G) 

(A.14) 

This form is equivalent to (jA.lip . but easier to work with 
as it does not require a sum over the set of paths on G. 

We point out that in this argument the sum over paths 
was only used to arrive at the form (jA.10[) for arbitrary 
A. Once this is known, the expressions (|A.14[) for the 
coefficients were obtained by Mobius inversion with no 
further reference to paths. This suggests that a shorter 
derivation may exist. 



2. Extension to n-point connectedness functions 



The expansion (|A.9[) generalizes readily to n-point con- 
nectivity functions; the criterion is simply that n root 
points xi , ■ • ■ , x„ are connected if and only if there ex- 
ists at least one path from Xi to each x;, 2 > ? > n, where 
we select xi arbitrarily. Note that in enumerating the set 
of paths from xi to Xi, we must include those paths that 
pass through other root points. Using inclusion-exclusion 
(|A.2p and equation (jA.6[) again, we may write the indi- 
cator function for this event as 



I[xi, • • ■ , x„ connected at p] — Ic(xi, Xi)|p 



1=2 



i=2 0cr;cr^i,>ci '-7er;. 



(A.15) 



Repeating the previous derivation and grouping together 
terms that test the same set of edges, we obtain the dia- 
grammatic expansion 

Cx,,...x„(p)= E d{G)Pr[G<p], (A.16) 

where the n-point d-weight is 
d(GeGx,,..,xJ = 

n 

n E (-l)l^*l+'l[U,r^ covers G]. (A.17) 

j=2 0crjcr^i,x,(G) 

Again, the fact that d{G) may be computed in terms of 
sets of paths covering G establishes that it is a topological 
invariant, unchanged by adding vertices of degree two to 
G. 



The argument following (jA.lip also carries though, 
since the above definition of the c?-weight may be ex- 
tended to arbitrary graphs and we may perform Mobius 
inversion on the connectedness function evaluated on an 
arbitrary n-point graph, obtaining 

d(GGGx,,..xJ= E (-1)1^(^)1-1^'! 

E'iZE{G) 

X I [E' connects xi, • ■ • x„] . (A. 18) 



Our final results, equation (|A.16p with (|A.18p . constitute 
a complete low-density expansion for all connectedness 
properties of percolation clusters. 



3. Equivalence with the Potts model 

The development of the field theory for the Potts model 
is described in detail elsewhere [13, [23, [1^, \^ and we will 
recall only the parts of the derivation that are relevant to 
our discussion here. The Q-state Potts model on a graph 
A [2l| has, associated with each vertex x of A, a degree 
of freedom q;(x) which may take on any of Q discrete 
states ( "colors" ) . The Hamiltonian for this model in the 
absence is 



H 



J E 



(x),q(x') 



1) 



(A.19) 



where the sum is over edges indexed by the two incident 
vertices x, x'. The partition function can be expanded 
in the form 20] 



E 



-f3H 



{q(x):xGI^(A)} 

= E n [(l-e"^'')'5a(x),a(x')+e 

{a} (x,x') 

= E pl^'l(i-p)i^i-l^'lQ'^=('= 



(A.20) 
-^•'] (A.21) 
(A.22) 



E'CE 



where p = 1 — e~^'^ and again E — E{A). When Q — > 1, 
the partition function becomes Z = 1, and the expan- 
sion corresponds to the sum of probabilities for the sets 
E' of occupied edges in bond percolation with indepen- 
dent probabilities p for occupying each edge. The Q-state 
Potts model partition function, viewed as a function of 
Q and p, is also (essentially) the Tutte polynomial [23| . 

The states at each vertex can be represented by an 
overcomplete set of Q vectors e*^, a = 1, . . . , Q, in a 
Q — 1 dimensional space. These vectors are obtained 
by projecting the position vectors of a regular Q-simplex 
in Q-dimensional space onto the subspace orthogonal to 
the vector (1,1, ...,1). More concretely, if we let the 
coordinates of these vectors with respect to some basis 
be ef , i = 1, . . . , Q — 1, the set of vectors may be uniquely 
defined up to relabeling and change of basis by requiring 
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that 

Q 

E < - 0' (A.23) 

Q 

J2efef = Q6,,, (A.24) 

Q = l 

Q-1 

and E efef = Q5a/3 - 1. (A.25) 
1=1 

In equations (|A.23P - (|A.25jl . we have normahzed the 
vectors following the convention used in [H, [2§, |4l| . Note 
that and [13] adopt a different normalization. 

To obtain the two-point connectedness function, we 
introduce factors e"j'^^\ ^'i^^^^ i^ito the sum. If they are 
not in the same connected component in the expansion, 
the sum over all a's gives zero by (|A.23p . That is, 

Cm,J2(x1,X2) = 

{a} (x,x') 

(A.26) 

= ^iiM E ^ connects Xi,X2] 

E'CE 

xpl^'l(l-p)l^|-|^'lQ^=(G-'), (A.27) 

where we also used (jA.24l) . After removing the factor 
(5i^.i2 and setting (3 = 1, this is equal to Cx,y(A,p). 
Now we rewrite 

P^Q(x),a(x') + (1 - P) = P(<5a(x),a(x') - 1) + 1. (A. 28) 

(Although this does not explicitly involve the ef s, this 
choice is motivated by the form of eq. (|A.25[) as Q ^ 1; 
note that there are many similar expressions that become 
equal to this for Q = 1.) We expand the Potts correlation 
function Ci^.i2(xi,X2) using this decomposition for each 
edge, and then once more for (5a(x).a(x') — 1 on each edge. 
This yields 

C,,,,,(xi,x2) =(5,,,,, E P^'^'^dQiGE'), (A.29) 
£;'C£;(A) 

where 

dqiGE') — E connects Xi,X2] 

E"<ZE' 

X (_l)l^'l-|iJ"lQJVc(B"). (A.30) 

Removing 5i^^i^ and setting Q = 1, we recover the ex- 
pressions eq. (jA.12p and (|A.14p . The derivation can be 
readily generalized, at least to the 3-point connected- 
ness function. Hence Essam's diagrammatic expansion 
is identical term-by-term with the low-density expansion 
of the Potts model in the Q ^ 1 limit. 



APPENDIX B: PROPERTIES OF MST PATHS 

In this Appendix we present proofs of properties 
obeyed by paths on the MST. These properties are used 
in Appendix [C] below to construct a diagrammatic ex- 
pansion for the MSF path vertex. 



1. MST paths as geodesies 

We first define a minimax path between two given ver- 
tices on the finite graph A: a (self-avoiding) path is a 
minimax path for the pair of (distinct) vertices x, y if 
among all paths from x to y it has the lowest value of 
the most costly edge (among all edges on the path). That 
is, it is a minimum (over the set of paths from x to y) 
of the maximum (over edges on the path) of the cost of 
the edge. We note immediately that in general there is 
more than one minimax path for the given vertices, even 
though they must all share the same most costly edge 
(we assume that no two edges have equal cost). We say 
that a path 7 is a geodesic if, for all vertices w, z lying 
on 7, the subset of 7 which connects w, z is a minimax 
path from w to z. A geodesic passing through x and y 
is necessarily a minimax path for x and y. Applying the 
definition of geodesic for the case where w, z are adjacent 
vertices connected by a single edge, we see that we may 
equivalently define geodesies as those paths all of whose 
edges are minimax paths connecting the vertices to which 
they are incident. We note that a geodesic cannot be a 
cycle, so it must have endpoints. We may also remark 
that the geodesic path is the correct strong disorder limit 
of the optimal path that which minimizes the total 
cost of all edges on the path with fixed endpoints. 

It is not always well-appreciated in the literature 
that a minimax path is not unique, see for example 
[lOl [H, [1^ , which frequently refer in the singular to 
"the" minimax path between two points. These sources 
really mean the geodesic path, which we will now prove 
is unique. 

We now prove that there is a unique geodesic between 
any two given vertices, say x, y (provided they are on 
the same connected component of the underlying graph) , 
provided that all edge costs are distinct. Specifically, we 
let 7 be a geodesic and we will show that no other path 
7 7^ 7 which shares the same endpoints may also be a 
geodesic. We noted above that any minimax path from 
x to y, such as 7, must pass through the same most costly 
edge ei. However, at this stage it is not clear that they 
all do so in the same direction. But if we consider the 
endpoint of ei that is encountered first on walking along 
7 from X to y, say w, then any minimax path, such as 
7, from x to that endpoint must pass through the same 
edge 62, which is the most costly on the subpath (but 
clearly less costly than ei). (If w = x, then we can start 
from y instead, and if Ci has endpoints x and y then 
we are done.) Note that this shows that 7 traverses the 
edge ei in the same direction as 7, because otherwise the 
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most costly edge on the minimax from x to w would be 
ei- Using induction on the number of steps on 7, we find 
that 7 must be the same as 7. 

By elementary properties of MSTs, all paths on the 
MST are geodesies. Likewise, all geodesies are contained 
in the MST, because in particular each of their edges is 
minimax for its two incident vertices, which is a property 
of the MST. Indeed, the MST of a graph is the union of 
all of its minimax edges. 

It is amusing to realize that the MST has the ultra- 
metric property (the content of this paragraph will not 
be used elsewhere in the paper). Let us assume that the 
costs are non- negative (if not, we can add a positive con- 
stant to all of them) . Then we can obtain a notion of dis- 
tance, or metric, d between any two vertices on the graph 
A, by defining i9(x, y) to be the largest (or minimax) cost 
on a minimax path from x to y, with 9(x, x) = if x = y. 
By definition, a metric should be finite and non-negative, 
symmetric (9(x, y) — 9(y,x)), equal to zero if and only 
if X = y, and obey the triangle inequality. The first three 
properties are clear, while it is easy to see that d obeys 
the stronger property that, for any x, y, z, 

9(x, y) < max (a(x, z), a(z, y)) . (B.l) 

These four properties imply that d is an ultrametric. 
Note that the ultrametric inequality eq. (jB.ip implies the 
triangle inequality. For ordinary metric spaces, one de- 
fines geodesies to be paths of shortest "length" using the 
metric, and this motivates our terminology above. Fur- 
ther, if X, y, z are three distinct points, the ultramet- 
ric property implies that if 9(x, y) < 9(x, z) and 9(y, z), 
then 9(x, z) = 9(y, z). It is well known that an ultramet- 
ric space with a finite number of points can be viewed as 
a tree, which we imagine depicted with the points as the 
leaves located on a hyperplane, other vertices to one side 
of the hyperplane, connected by straight lines, and the 
ultrametric represented by the height above (in the di- 
rection orthogonal to the hyperplane) the leaves to which 
one must go in walking from one leaf to another along 
the tree. In the present case, this essentially corresponds 
to the MST. The tree is trivalent (except at the leaves) 
with probability one. The trivalent vertices represent the 
edges on the MST, with their height as their cost. In fact, 
if we consider the subforest of the tree consisting of the 
vertices at height less than or equal to some bound, then 
this represents the MSF(p). 

2. Identifying MST paths through binary 
comparisons 

Let 7mst(x, y) be the geodesic from x to y, or equiv- 
alently the path on the MST. The geodesic property of 
7mst(x, y) allows it to be selected from the set rx,y of all 
paths connecting x, y by means of repeated comparisons 
using a binary ordering relation defined as follows. 
Let 7 and 7' be two paths in Fx^y. Let ei, 62, • ■ • e„ be 
the first, second, ... n-th most expensive edges on 7, and 



likewise for e'^, 62, ■ • ■ on 7'. We say 7 ^ 7' if and only 
if there exists some j such that < ^e'. and i^i = ^e'. for 
all i < j: in other words, we compare the most expensive 
edges whose costs are not identical. We will prove that 

7MST(x,y)== niin 7; (B.2) 

in other words, 7mst(x, y) is the minimal element of the 
set Fx,y under the ordering defined by -<. 

Again, we assume all edge costs to be distinct, which 
implies that = e[ <J=> = f-e'. and -1(7 < 7') <S=^ 
7 )^ 7': i.e., the relation -< defines a total order on the 
set of all paths between fixed endpoints. Under this as- 
sumption, for any two paths 7, 7' G Fx.y we have 

7^7' max le < max 1^' ■ (B.3) 

eG7 — 7n7' e' ^-y' —^r\^' 

Let 7 be a path satisfying 

7 = min 7'. (B.4) 

In particular, 7 is less than all paths in F(7; x, y) C Fx,y, 
the set of paths in Fx,y having no edges in common with 
7. By (jB.3p . this means -< compares only the most ex- 
pensive edges on the paths, so if 7 -< 7' G F(7; x, y) then 
7 must be a minimax path for x and y. 

Similarly, for any subpath 7 C 7 with endpoints w, z, 
we define F(7; w, z) C Fx,y as the set of those paths that 
coincide with 7 from x to w and from y to z, and have no 
edges in common with 7; in other words, the set of those 
7' such that 7—707' = 7. Because 7 ^ 7' G F(7; w, z), 7 
is a minimax path from w to z. This holds for all choices 
of w and z, and hence 7 is a geodesic. It is unique and 
can be identified as a path 7 = 7mst(x, y) on the MST 
by the results above. 

We finally note that identifying MST paths through 
the definition (jB.4p is most convenient for the purposes 
of our diagrammatic expansion (jII.18|) . It would be very 
inefficient computationally, since we make many unnec- 
essary comparisons with paths that are not in any of the 
sets F(7;w, z). On the other hand, the geodesic char- 
acterization of MST paths is less directly useful for our 
purposes, since it requires keeping track of the locations 
of the most expensive edges. The geodesic properties of 
MST paths are very useful computationally: they are es- 
sential in constructions of linear-time algorithms for MST 
path verification lj:6] which were used in |47j] to give a 
randomized algorithm which constructs the entire MST 
in linear time. 



APPENDIX C: DERIVATION OF LOW-DENSITY 
EXPANSION FOR MSF PATHS 

In this Appendix we derive the exact low-density ex- 
pansion for the probability C^y{p) that the points x, y 
are connected by a path on the MSF(p) which passes 
through z, on a finite graph A. 
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First, we can formally define an indicator function 



Ic(x,y;z)lp 



and then 



[(x, y connected at p)A 

(7MST(x,y) passes through z)] (C.l) 



C^,y(p) = (I,(x,y;z)|p). 



(C.2) 



In the Kruskal process, edges are never removed from 
MSF(p) as p is increased, so if a path connecting two 
points on the MSF at a parameter value p exists, it must 
be identical to the unique path connecting those points 
on the completed MST. In Appendix IB 21 above we use 
the geodesic properties of the MST path to arrive at the 
definition 



7mst(x, y) = min 7, 



(C.3) 



where min denotes the minimal element under the rela- 
tion -< defined in (|B.3|1 . This lets us write the indicator 
function in (jC.ip as a sum over all paths 7 G ^x,y, in the 
form 



Ic(x,y;z)|p = 



E 



I[7 < p]4^)[7]I 



A (7^7') 



where 



I^^^ [7] = I[7 passes through z] 



(C.4) 



(C.5) 



where again Fx^y is the set of all paths on A with end- 
points X and y. For the time being, we will suppress the 
dependence of all expressions on the underlying graph A. 

Expanding the indicator function I A-f'er — 7') 
by inclusion-exclusion gives 

Ie(x,y;z)|p- J2 lb<p]t^b] 



^ (-l)ir'll /\ -(7^7') 



(C.6) 



Because the uniqueness of the edge costs implies -1(7 ^ 
7') <;=^ (7 >- 7'), we may restrict the sum over subsets 
of Fx,y to those not containing 7 itself. We now reorga- 
nize the double sum by grouping together all terms that 
test the same set of edges, as was done for equation (jA.9[) 
for percolation. For each term in (|C.6p . the edges in 7Ur' 
form a graph G in the set Gx.y of all vertex-irreducible 
graphs with root vertices x, y. When we regroup the sum 
in terms of these graphs, we obtain a sum over sets F' 
of paths from x to y which cover G (these sets are the 
previous 7UF' redefined as F', so contain the chosen path 
7), similar to what was obtained to percolation. Unlike 
the percolation case, we still have the outermost sum in 



(|C.6|) . which becomes the innermost sum over elements 
7 of F'. Thus the expansion becomes 



Ic(x,y;z)|p = 



GeG^ 



E 

r'cr^,y(G) 



E 

7er' 



(_l)|r'|+i][[r' covers G] 

tHl]- (C.7) 



A (7^7') 

7'er' — f 



As in the derivation of (jA.lOp . we may factor out the 
dependence on the parameter p as 



Ic(x,y;z)|p = 



E 



^MSF 



(G)I[G<p], (C.8) 



where we have introduced dMSF(G), the analogue of Es- 
sam's c?-weig ht (|AjT|i for MSF paths: 



duspiG) 



J2 (-1)I^'I+'I[F' covers G] 
r'cr^.y(G) 



xE^ A (^^^') 

7er' L^'er' — 1 



4==) [7]. (C.9) 



'^msf(G) depends implicitly on x, y, and z, and on the 
costs of the edges of G. 

In the analogous statement (jA.lip for percolation, we 
found d{G) was independent of edge costs. Here to eval- 
uate dyispiG) we need to compare paths which cover G 
using the relation From the definition (|B.3p . a nec- 
essary and sufficient set of information to do this is the 
relative ordering of the edge costs of G. As discussed in 
Section III Bl we will introduce an ordering by indexing 
the set E of edges of A arbitrarily and defining an order- 
ing of their costs to be given by a permutation tt G S\e\ 
on the set of \E\ elements, via 



I < J. 



(C.IO) 



Then the induced ordering on a subset E' of E is written 
tte'- With this notation, we see that dMSF is a function 
of the graph G and edge cost ordering tte{g)^ so we write 
rfMSF(G|7r^(G)) (it still depends implicitly on x, y, and 
z). 

The second sum in (jC.9l) detects whether the maximal 
path in F' passes through the point z, so for a fixed edge 
cost ordering tt we may write 

dMSF{G\TrE(G)) = 

J2 (-1)I^'I+'I[F' covers G]4")[max7]. (C.ll) 

r'cr^.y(G) 

Note that, as a consequence of our use of inclusion- 
exclusion, this result is mildly counterintuitive: we are 
attempting to calculate the probability that the MSF 
path passes through z, and by definition the MSF path 
(if it exists at p) is the minimum out of all paths in Fx.y. 
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However, for each graph in the expansion of this proba- 
bihty, the relevant event is that the maximal path of the 
covering F' passes through z. 

We may now take the expectation value of Ic(x, y; z)|p 
over all realizations of the edge costs in order to obtain 
the analogue of (|A.10|) . 

^x,y(p) dMSF{G\-KE{G)) 

xPr[7r£;(G) A(G<p)]. (C.12) 

We note that the ordering t^e(g) ^-nd the event that all 
edges of G be less than p are independent, so the last 
probability factorizes. Because the edge costs are iid, 
all orderings of the edge costs are equally probable and 
Pr[7r] = l/\E{G)\\. 

We may find an alternative expression for cZmsf in 
terms of a sum over edge subsets instead of sets of 
covering paths, analogous to our derivation of (|A.14p 
from (jA.ll[) . The argument proceeds the same way: we 
first make explicit the dependence of MSF path con- 
nectedness functions on the graph A in the set Gx,y,z 
of graphs containing the three root points x, y, and 
z, writing it Cxy(A,p). Because (jC.lip contains a 
factor of I[r' covers G] in the summand, we also have 
'^MSF(G|7r) — for disconnected or vertex-reducible 
graphs. The sum in (jC.12|) may therefore be extended 
to all subgraphs of G as 

TrGS|B(A)| E'CE(K) 

xdMsAGE'WE'mGE' <p]Pr[7r£;.]. (C.13) 

Because dMSF(G\i: e(g)) is dependent on tt, we must work 
under the sum over edge cost orderings in performing 
the Mobius inversion step. We therefore work with the 
conditional quantity 

G^,y{G,p\T^E(G)) = E '^MSF(G£;'|7r£;/)Pr[G'£;' < p] 
E'(1E{G) 

(C.14) 

appearing as a summand in (jC.lSp . Evaluating this at 
p = 1 yields 

l[E{G) connects x, y]4^) [7MST(G|7r£;(G))] = 

E '^MSF(G_E'|7r£;'), (C.15) 
E'<ZE 

where 7MST(G|7r) is the path connecting the root points 
X, y on the minimum spanning tree of G obtained under 
the edge cost ordering tt. Mobius inversion of this sum 
gives 

E'QE(G) 

X l[E' connects x, y]4") [7mst(G£;' I^bOI- (C.16) 



This definition of c^msf is more convenient than (jC.lip 
for the proofs of appendix ID 41 It is also, in principle, 
more convenient for computation, since for large graphs 
the size of the set rx.y(G) of self-avoiding walks grows 
faster than \E{G)\, hence the sum in (jC.16p is more easily 
performed than that in (jClip . 

APPENDIX D: RENORMALIZABILITY OF THE 
MSF PERTURBATION EXPANSION 

In this section we give the proofs outlined in Section 
IIIIE| which establish that our perturbation expansion 
for MSF paths is renormalizable. Recall that diagrams 
of this theory with no path vertex are identical to those 
of percolation theory and hence pose no problem, while 
we construct diagrams involving the MSF path vertex by 
the substitution (jIII.7P : 

d{g)i{g) ^ 

E rfMSF(ek)OMSF(7r,to)/(e,{t4), (D.l) 

•"■eSiEcs)! 

where the integrals I{G) and I{Q, {te}) are identical Feyn- 
man integrals with only cubic interaction vertices, con- 
taining the factor for each such interaction, but in the 
latter integral the mass-squared is generalized to a dis- 
tinct parameter for each edge e of the graph Q. In this 
Appendix, we will drop the prime from the orderings tt' 
throughout; orderings T^E(g) are nonetheless the induced 
orderings on the set of highest costs of the set E{Q) 
of edges e € E{Q) oi & topological graph (Feynman dia- 
gram) or subgraph Q. dyisF and Omsf were defined in 
(111.81) , piI.Sp respectively, and as defined both depend on 
the structure of the entire graph Q. In particular, it is not 
entirely evident from the definition (jll.Sp how dMSFl^l'"') 
could be computed from knowledge of its values on sub- 
graphs of Q. 

The Appendix is structured to give proofs of the fol- 
lowing results. We begin in section ID II by introducing 
terminology common to all sections of this Appendix and 
explaining the parametric formulation of Feynman inte- 
grals. In section ID 21 we obtain the effect of the Omsf 
operator, which is very simple in the parametric formu- 
lation: it introduces a simple product factor i^^ into the 
integrand, which depends on the choice of an ordering tt 
for the costs on the graph. 

In section lTTSl we prove that the superficially-divergent 
subintegrations (as the cutoff A oo) associated to a 
connected subgraph Ti. come only from a subset S' C 
S\E(g)\ of all possible orderings on the edges of Q. Specifi- 
cally, if 7"^ is a three-point subgraph or a 2-point subgraph 
containing the path vertex, then 5" = ttjo] , in which all 
costs in the subgraph are cheaper than all those out- 
side. Similarly, if 7i is a two-point subgraph (i.e. a self- 
energy) , the only superficial divergences are for orderings 
S' = 7r[o] U 7r[i] , in which at most one edge in Ti, has cost 
higher than one or more outside Tl. Moreover, for these 
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orderings, with one class of exceptions the divergences in 
self-energy or cubic coupling (3-point) subgraphs are the 
same as those in the corresponding percolation diagrams, 
up to the d-weights. These results generalize easily to di- 
agrams with several superficially-divergent subdiagrams, 
if these are pairwise either disjoint or one inside another. 

Having identified the important orderings, we consider 
in section ID 41 the behavior of the dMSF weights for these 
orderings. We show that the weights obey nice factoriza- 
tion properties for connected subdiagrams with two or 
three external points for orderings in class tt^qj , and also 
(in a different, more general form) for self-energy sub- 
diagrams H with orderings in which one or more edges 
in Ti is more costly than at least one outside Ti. The 
factorization has the form 

dMSF{Q\T^E(G)) ^ d{7i)dMSF{Q l'H\T:E(g /H)) (D.2) 

if the path vertex is not in TL (note the appearance of a 
d- weight from percolation), and 

dMS,F{Q\T^E{g)) = dMS,F{'H\-K E(n))dMSF{Q l'H\-K E(g /H)) 

... 

if the path vertex is in Ti. The precise definitions, in 
particular for T^E(g/H)) will be given in section fP 41 Here 
and below we use notation Q/Ti to denote the diagram 
obtained by contracting the subgraph 7i to a single ver- 
tex (which may be of degree 2, producing a harmless 
extension of the class of diagrams to be considered) . 

In section [D3] we come to the heart of the proof. We 
use a theorem of Bergere and Lam [s^l to show that the 
Feynman integral for each diagram in our perturbation 
expansion can be rendered absolutely convergent by a 
procedure of subtracting all the superficially divergent 
parts of the integrand. Furthermore, utilizing the results 
of preceding sections, all the terms that have to be sub- 
tracted for divergent subdiagrams (including those con- 
taining the path vertex) are the same as those for a cor- 
responding full diagram, with the exception of one class 
of terms as mentioned above, which is dealt with in sec- 
tion ID 61 Subject to the latter result, this means that 
all divergences are dealt with by renormalizing parame- 
ters and the overall scale of the vertex functions, as in a 
renormalizable field theory. 

Finally, in section ID 61 we prove that the class of ex- 
ceptional subleading divergences in the self-energy sub- 
diagrams cancel in the sum over all diagrams of a given 
order. This completes the proof of renormalizability to 
all orders in the perturbation expansion. 

1. Definitions 

We begin by considering an arbitrary Feynman integral 
associated with a diagram Q appearing in the perturba- 
tive expansion of a correlation function in, for example, 
the field theory of percolation. 

First, we recall the expression for a Feynman integral 
associated with Q given in eq. (|HI.2p . Let V, E be the ver- 




FIG. 3: Depiction of an arbitrary graph contributing to the 
MSF path connectedness function. The various root vertices 
used in the definition of this function are labeled with open 
circles. 



tex and edge sets of and let N'e.v be its incidence matrix 
under an arbitrary orientation of its internal edges; i.e. 

A/'e,« = 1 if w is the head of e, 

= — 1 if ti is the tail of e, and (E).4) 
= otherwise. 

Neglecting the cut off for a moment, the integral with 
which we are concerned is 




(D.5) 



Here i}^cy±)v is the net external momentum incident on 
the vertex v. Since we will replace the percolation d- 
weight with the appropriate dMSF-weight, we neglect the 
factor d{Q) and also to the power of the number of 
internal cubic vertices (or other couplings for interactions 
of different degree that may be present more generally). 

We make progress by expressing the Feynman integral 
(ID.Sp in terms of integrals over the Schwinger parameters 
a. This makes use of the identity 1/X = J^da e""^ to 
rewrite part of the integrand as 

(For brevity, we let Ag denote the set of parameters 
introduced above.) The total of d\V\ (5-functions can be 
rewritten using the identity 2tt 5{k) — J d\e^^^ for each. 
Integrals over the internal momenta are now Gaus- 
sian and can be performed, and then the A-integrations 
become Gaussian and can be performed, except for one 
which produces a i5-function expressing conservation of 
the total momentum, (27r)''(5''(^^gy (kcxt))j)- Omitting 
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this 5-function, we have 



1 

X exp -k^^tA^^koxt - ^ a^t^ . 



(D.7) 



Here A is a \V\ x \V\ matrix, which is a Laplacian on 
defined by 

^ — ^ rv . 



and i'(^e) is defined as 



P{Ag) - n "0 det'A(Ag) 



(D.9) 



\eeE 



in which the determinant det' is that of A with one row 
and column removed, so as to remove the zero mode, 
kcxt is viewed as a -component vector, and C = C{Q) 
is the cyclomatic number of the number of independent 
loops (cycles) of Q. P{Ag) is a homogeneous polynomial 
of degree C. These expressions are quite general and may 
be obtained for the diagrams of any field theory; for a 
further discussion consult [13, ■ It is interesting that 
P{Ag) and also A~^ can be related to weighted sums 
over spanning trees on Q [H, [s^ by the Kirchoff matrix- 
tree theorem [3]; it is not clear to us whether this fact 
is deeply involved in the renormalizability of the theory 
of MSF(p). 

For many diagrams, the integral I{Q) as written in 
(ID.SP or (|D.7p is ultraviolet divergent and must be regu- 
larized, which is done by restricting the momentum inte- 
grations in (|D.5|1 to the region |qc| < A. We implement 
this in (|D.7P by taking the range of integration of each 
of the to be [A~^, oo), which exponentially suppresses 
contributions from jq^l ^ A. 

We make use of the parametric representation for Feyn- 
man integrals for two reasons. First, it greatly simplifies 
the study of renormalization of the expansion, as in the 
field theories in [H, (S^, 113 ■ 

The second reason we employ the parametric repre- 
sentation is that, as we show in the next Subsection, the 
action of Omsf takes a particularly simple form. Ap- 
plication of Omsf directly to (jD.5[) results in intractable 
integrals over the {t^} for diagrams beyond one-loop or- 
der, while we are able to obtain its action on an arbitrary 
graph in closed form in equation (jP.lip . 



2. Effect of Omsf operator 

The preceding discussion applied to the diagrams from 
the field theory for percolation. To investigate how things 
change when we calculate MSF diagrams, we specify a 
given total ordering tt of the masses of Q, such that 



i < j 



(D.IO) 



The diagrammatic contribution to the MST theory is ob- 
tained by summing over all total orderings of edge costs 
consistent with the placement of the path vertex, accord- 
ing to (|C.16p . To find the contribution from one ordering 
TT, we apply the operator C'MSF(7'',t) defined in (jIII.5[) to 
both sides of (|D.6|) . obtaining 



00>t^(l) >••• 
■■■>t„(lEl)>t 



eeE 



d 1 
dt^ q2 -I- 



(D.ll) 



°° -pr «7r(i) C^"7r(i) ajq;+f) 
^=l Ej = l 



The integrand on the right-hand side is that appearing 
on the right-hand side of (|D.6p . multiplied by a factor 



\E\ 

n 



= 1 



(D.12) 



(We note that FTr{Ag) is of the same form as we obtained 
on the lattice in equation (jll.isp . This is another man- 
ifestation of the well-known equivalence between scalar 
field theory and a system of random walkers.) Thus fi- 
nally our prescription for evaluating the contribution of 
each diagram is that it is given by the parametric Feyn- 
man integral as for percolation, but with the factor 



J2dMSF{g\n)F^{Ag) 



(D.13) 



inserted inside the a integrals, replacing the d(Q) weight 
for the percolation theory. 

The factor F^{Ag) obeys < F^(Ag) < 1 for any 
Ag e [0, oo)l^l, and has the property that it reduces to 
one as we go towards the limit in which 



<C a 



TV{\E\ 



(D.14) 



It tends to suppress orderings which do not obey the ver- 
sion of these inequalities in which all ^'s are replaced by 
<'s. Thus it acts to replace the strict inequalities on the 
te's by corresponding but softer conditions on the a^'s. 
This result makes intuitive sense: high-momentum (small 
a) propagators correspond to lattice walks consisting of 
relatively few edges. In the Kruskal process, we expect 
the shortest paths to be completed first, at the lowest 
value of p, corresponding to a larger mass-squared tQ. 

It will be useful to simplify the factors as much 
as possible, by performing (or partially performing) the 
sums of dMSF{G\TT)F.^ over orderings tt as much as pos- 
sible before performing the integrals. We now give some 
basic formulas that are a step in this direction. First, we 
obtain another proof of (|II.17p from the fact that 



(D.15) 



A more general fact that will be useful is that if we con- 
sider a subset of edges E' <Z E oi E and orderings tt 
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such that the masses on edges of E' are greater than all 
those in E — E' , and sum over all such orderings that 
fix an ordering on E — E' (such orderings can be writ- 
ten as TT = CT o TTo for ttq any one such ordering and a a 
permutation in Se> C Se), then: 

o-GSe/ i=\E'\ + l ^j=l^^o{j) 

in which the right-hand side is independent of the choice 
of ttq. This follows by using eq. (|D.15p applied to the 
restricted sum over orderings. Indeed, as the derivation 
of this identity only used the sum over a smaller set, this 
can be used in a proof by induction (on the size \E\) of eq. 
(|D?T5l) itself. The induction step, of ta king \E '\ ^\E\~1 
and summing the right-hand side eq. (|D.16p over cosets 
5'|_E|/'S'|E|_i is simple. 

3. Estimating MSF Feynman integrals 

In this Subsection, we describe how the divergent be- 
havior of a given diagram of the MSF path theory differs 
from that of the diagram from percolation theory from 
which it was obtained, and obtain some basic statements 
about the form of the divergences for each ordering. 

In the absence of the F^^ factor, the parametric form of 
the Feynman integrals may in general suffer from diver- 
gences associated with the region a — > for some or all 
a's. These take the place of the possibly more familiar 
divergences at large k in the original momentum space 
integrals over q^; recall that the latter integrals have al- 
ready been done, after exchanging orders of integration. 
For a IPI graph Q, the superficial degree of divergence of 
I{G) is obtained easily from the momentum-space form 
by counting the total number of powers of all q^'s and 
integrations J d'^q^, and is given by 

uiG) ^ dCiG) - 2\Eig)\. (D.17) 

This formula holds for any field theory of scalar fields in- 
teracting via non-derivative couplings. The same result 
is easily obtained in the parametric representation also 
[33I [35I [sgI \3T\ . It may be obtained more formally by 
rescaling the p^a^ for all edges of Q, with p ^ 0+. 

The formula may also be applied to the subintegral asso- 
ciated with a connected IPI subdiagram TL oiQ (strictly, 
a subdiagram is a subset of the vertices of together 
with all edges that connect these vertices); this will be 
denoted uj{H). In this case, it is obtained from the be- 
havior as the subset of ae associated with edges of H are 
scaled to zero by a common factor. Notice that the su- 
perficial degree of divergence for a subgraph Ti. might be 
larger than that for Q. A graph or subgraph is said to 
be superficially divergent if its superficial degree of diver- 
gence is positive or zero, and superficially convergent if 
its superficial degree of divergence is negative. (A graph 
with Lu = may diverge more slowly than any power of 



A, for example logarithmically, or may be convergent.) 
It is a theorem that if oj is negative for Q and for all its 
subgraphs, then the associated Feynman (parametric) in- 
tegral is absolutely convergent. 

For the theory with cubic interactions that we consider 
here, the only connected IPI graphs that are superficially 
divergent at d = 6 dimensions are (a) any self-energy 
diagram (with two external points), because all have lu = 
2, (b) any vertex correction diagram, that is a graph with 
three external points, because all have w = 0, and c) a 
self-energy graph with a (jx^ insertion, which have the 
same form as the vertex diagrams in b). Here in a) and 
b) an external point means that a line that "leaves" the 
graph (joined to it by a cubic vertex like the others) was 
removed to leave the IPI part. The graphs containing cf)'^ 
are relevant to the path vertex that we wish to consider 
in this paper. Other graphs are superficially convergent. 

Turning to our theory for MSF(p), the parametric form 
of the Feynman integral for a given ordering tt is simply 
modified by the insertion of the factor ^^(Ag). Because 
FTr{Ag) is bounded, it follows that the superficial diver- 
gence of any diagram or subdiagram of the MSF theory is 
no worse than the corresponding diagram of percolation 
theory from which it was obtained. More formally, F-^ is 
a homogeneous rational function of degree zero, and so 
the superficial degree of divergence for Q is again 

However, for a subgraph H of Q, F^(Ag) may reduce 
the superficial degree of divergence below uj{TL). RecaU 
that for a subdiagram, we consider the limit as for 
e G E{Ti.) go to zero simultaneously, by scaling them with 
a common factor p^, leaving for e ^ Q — Ti. unchanged. 
Considering each of the |-E(W)| factors in F-^{Ag) that 
have numerator for an edge e G E{Ti.) appearing in 
(|D.12p in this limit, we see that in this limit p ^ 0, 

F,(Ae) = o(p2"'(«'^)), (D.18) 

where we define nT^['H, Q) to be the number of masses 
for e G E{Ti) that are less than at least one of the masses 
in E{Q) under the ordering tt. (Clearly n-^{T-L,Q) = for 
Ti = Q.) For a fixed IPI connected subgraph Ti, of Q, this 
provides a useful partitioning of orderings into sets 7r[,„] , 
m — 0, 1, . . . : 

■K[m\ = {tt : n-^iU, Q) = m}. (D.19) 

Thus the orderings 7r[o] (which will prove most important 
in what follows), for which = 0{1) as p ^ 0, are those 
where all of the masses on the edges of the subgraph Ti. 
are larger than those in C/ — Ti, that is all the costs in TL 
are lower. 

We may add this result to the superficial degree of di- 
vergence to obtain the overall superficial degree of diver- 
gence of a connected IPI subdiagram of a connected 
IPI diagram Q under the ordering tt: 

wmsf(h, g\Ti) = uj{n) - 2n^{n, g). (D.20) 

This implies that it is only for class ttjq] that the superfi- 
cial degree of divergence of the subgraph Ti. is unchanged 
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by F^. For subgraphs with Lu{Ti,) = (i.e. the vertex 
and path vertex diagrams), orderings other than those 
in 7r[o] give convergent subintegrals. For the self-energy 
subgraphs, with to — 2, orderings in ttjx] lower the super- 
ficial degree of divergence to 0, and these are additional 
divergences with which we will have to deal. Moreover, 
in all cases there are subleading terms in the behavior of 
as p ^ for a subgraph, and while these terms are 
superficially convergent in most cases, the first sublead- 
ing term also has zero superficial degree of divergence in 
the case of the self-energy subdiagrams. 

For further analysis, it is helpful to consider the sum 
X^TT '^MSF(0|7r)F7r(^g) and to attempt to simplify it as 
much as possible, so that the evaluation of the paramet- 
ric integrals reduces to those for percolation as much as 
possible. Indeed, by the "contribution of a diagram" in 
general we mean the weighted sum over orderings. In or- 
der to consider divergent subintegrals for subdiagrams, 
it is useful to have factorization properties of the weights 
c^MSF- It is to this that we turn next. 

4. Factorization properties of MSF diagrammatic 
weights 

In this section we demonstrate that the diagrammatic 
weights diviSF possess enough factorization properties for 
our proof of the renormalizability of the perturbation ex- 
pansion to go through. Let us first recall that for the 
d-weights in percolation, the weight for a diagram Q con- 
taining a 2- or 3-point subdiagram Ti factors into the 
weight for H times that for the "quotient graph" Q/H 
in which the subgraph Ti is shrunk to a single vertex 
(formally, its vertices are identified, and its edges are 
deleted): d{Q) = d{T-L)d{Q /Ti). This is immediate in the 
Potts model formulation in which the d- weights originate 
from contracting together tensors, due to Sq permuta- 
tion symmetry (apart from the problem of giving a formal 
definition of the Q ^ 1 limit). It can also be derived from 
the combinatorial definitions described in Section [X] (this 
is shown in the case of some 2-point subdiagrams in Ref. 
[2^). It is important for the proof of renormalizability, 
as the contributions of such subgraphs in Feynman inte- 
grals will be treated as "correcting" or "renormalizing" 
the parameters attached to 2- and 3-point vertices in the 
graphical expansion. We require some similar properties 
in the expansion for MSF(pc)- 

Recall that the weights can be defined as in eq. (jC.16p 
[iovE^E{g)], 

duSFiGln) = J2 (-I)"'!""''! 

E'CE 

X l[E' connects x, y]4^' [7MST(e^iS' kfi')]- (D-21) 

This differs from the diagrammatic weight for percolation 
(IA.14P only in the presence of the additional indicator 

function I'f'' [jmst{Ge'\t^e')]- A graph that is not IPI 
can be decomposed into (connected) IPI subdiagrams 



lying on a chain of single edges and such IPI parts that 
form a path from x to y, and possibly other IPI parts. 
That is, Q may be constructed as a tree Qo (with x, y, 
z marked) which is then decorated by replacing its ver- 
tices V with subgraphs Hy As the MSF path must pass 
through a chain of IPI parts, it follows that for those 
IPI subdiagrams that do not contain the vertex at z, the 

fz) 

indicator Ic is independent of the path through such a 
IPI subdiagram, and accordingly the duSF weight fac- 
tors into a product of weights for the single edges and for 
the IPI parts. Moreover the diviSF factor for each such 
IPI subdiagram reduces to d in percolation for that sub- 
diagram (for a single edge, the weight is 1). Likewise, 
for a vertex-reducible subdiagram or "tadpole" , such as 
a IPI part connected to the rest by a single edge, the 
^MSF-weight is the same as in percolation and vanishes. 
Similarly the c^msf weight for a diagram that contains a 
disconnected subdiagram vanishes. Hence from here on 
we need consider only connected, vertex-irreducible IPI 
diagrams G that contain the path vertex at z, as well as 
root points that we can relabel as x, y. 

For MSF(p), the weights c^msf depend on the ordering 
TT of the costs of the edges of the topological graph Q, as 
well as on C/. In this section, we will denote these costs by 
the original symbol £e for edge e € Q (these costs in fact 
stand for the maximum, earlier denoted L^, of the chain 
of edges that are the image of e under an embedding of 
Q in the lattice). In terms of the costs, the ordering tt is 
defined by 

i<j ^ 4w <e^U)- (D.22) 

(We use the costs, rather than the mass-squared's for 
which the inequalities are reversed, because the authors 
find that this aids their intuition about MSTs.) In seek- 
ing a factorization similar to that for the d- weights in per- 
colation, there are two issues. Because the divisF-weights 
depend on a choice of ordering, one issue is whether some 
factorization holds at all for each ordering, and a second 
is, if there is some factorization, what ordering would be 
used for the quotient G/H. What we obtain below may 
not be the most general possible result. Instead we ob- 
tain statements for two (overlapping) sets of conditions, 
and these are sufficient for our purposes. 

Motivated by the considerations of which orderings 
produce ultraviolet-divergent Feynman integrals associ- 
ated with a subgraph, we first show that for orderings 
TT such that all edges in a connected subgraph Ti have 
lower cost than all others in Q, where 7i is a 2- or 3-point 
subgraph, and in the 2-point case the vertex z can also 
be present, factorization holds: 

dMSF(^/k) = dusFi'H\'^n)dMSFiG /HlTTg^T-c)- (D.23) 

Here the right-hand side involves the ordering tt-^^, which 
is TT restricted to 7i, and irg-n which is tt restricted to 
Q — H. (For graphs G, subgraphs H, and quotients GfH, 
we will allow abuses of notation like t^ech) = "''H-) Fur- 
ther, in the case in which Ti does not contain z, we al- 
ready know that dMSF{'H\ni-L) = d{'H). We recall that 
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these orderings are those in class tt^o] , which produce the 
leading divergence for the 2- or 3-point subgraphs. Fur- 
ther, the factorization generahzes to the case in which 
there are several disjoint such subgraphs, and the costs 
in the union of the sets of edges of the subgraphs are lower 
than those in the remainder of Q (regardless of the rela- 
tive orderings among the edges in the subgraphs). In this 
case, each disjoint subgraph carries a weight as for the 
single subgraph considered above. Then, because the 2- 
point (or self-energy) subgraph (that does not contain z) 
also has subleading divergences that occur when its costs 
do not obey the preceding conditions, we also derive a 
more general result for such a subgraph for any ordering. 
These results can be combined to handle a large class of 
orderings and subgraph structures. 

First we show that, if 7i is a 2- or 3-point subgraph, 
then in the sum over subsets of edges E' C_ E in dusF we 
can replace 

I[E' connects x, y] = I[E'{H) connects {x^}] 

X I[E'{g/n) connects x,y], (D.24) 

where x, y are the root points of Q, and {xi} are the 
root points of H, because other terms cancel. To see 
this, first notice that if for a 2- or 3-point subgraph TC, 
the "diluted" edge set of H, E'{H) = E' n E{H) [and 
similarly for E'{g/'H)], does not connect all the root ver- 
tices, then there is at least one root vertex not connected 
to any of them (this does not hold for a subgraph with 
more than three root points). Choose one of these, and 
without loss of generality suppose it is xi. In Q there is 
a single edge incident on Xi that is not in Ti (call it e^^ ) . 
The minimum spanning tree path ^mst{Ge'\t^E') from x 
to y on E' clearly cannot pass though Xi for such an E', 
whatever the ordering tt. We can pair off such subsets E' 
by choosing pairs of E' which are the same subsets except 
that the edge is in one and not in the other. These 
subsets differ in size by one, and the indicator function 

Ic^'' [7mst(^_e' \t^e' )] takes the same value for both. Hence 
these contributions cancel, and the result follows. 

Now we turn to the factoring of I'f^ [7mst(^?_e' {t^e' )] ; it 
is here that the form of the ordering enters. The case in 
which all edges in the subgraph H have costs lower than 
all those in G—Ti- is quite simple. First, the same property 
is inherited in the ordering tte' restricted to E'. As the 
Kruskal process runs on E' , these edges are tested first, 
and when that is completed the root points of the sub- 
graph are connected (this follows because we have shown 
that E' connects these vertices). For the remainder of the 
process, from which the path 7mst(S£;' Kb') is obtained, 
the subgraph Ti [or its diluted version which we denote 
'He'(H)] can be viewed as collapsed to a single vertex to 
produce Qe' l'HE'{n)- It is useful now to distinguish two 
cases: either z is in 7i, or it is not. In the first case, 
the MSF path must enter TL to reach z, and then leave. 
This implies that a) on Qe' /'He'(t-i)j the image of 7i is 
the point through which the MSF path must pass, and 
b) once within Ti. the path must pass through z. That is. 



we can write for the indicator function 

4^^[7MST(t/£;'k£;')] = ^'c'\lMST{'HE'{-H)\T^E'{n))] 

X l^/'HlMSTiGE'/nE'iH)\T^E'-E'(n))] (D.25) 

The summation over subsets of the edges E' can be 
written as a sum over subsets E'{7i) and over E" = 
E' — E'{'H), and so the factorization of the dMSF-weights 
as in eq. (|D.23|) follows. Likewise, in the case where the 
path vertex is located in Q ^ Ti, wc can simply write 

I^c''bMST{gE'\7TE')] 

= li^^bMST{gE'/'HE'(H)kE'~E'(H))], (D.26) 

and again the form in eq. (|D?23l) follows, though now 
dMSFiTi-liTT-c) = d{TL). Together these prove all the rela- 
tions shown in Fig. 2] for the the stated class of orderings. 




FIG. 4: Factorization properties of dMSF(5|7r) for connected 
subgraphs that are proved in the text for orderings in the class 
vr[ol. 

For graphs with several disjoint subgraphs of the same 
type as the single connected subgraphs considered above, 
the same proof goes through if all edges in all the sub- 
graphs are less costly than those outside. That is, the 
edges in one of the connected subgraphs need not be all 
more, nor all less, costly than those in another one of the 
subgraphs. Thus this result is more general than simply 
iterating the application of the preceding result, though 
the final factored form of the dMSF-weight is the same as 
if it were. 

We now turn to a more general argument for the case 
of a connected IPI 2-point subgraph Ti, and any ordering 
TT. It holds if the path vertex z is within 7i, however for 
orderings not in the class 7r[o] which are already covered 
by the preceding proof, the corresponding Feynman in- 
tegrals are convergent, so we will not make of this, and 
can assume that z is not in TL. 

Let e-j-i be the most costly edge in TL for the ordering 
TT. Now we consider the evaluation of the dMSF-weight. 
For each diluted edge-set E' ^ this involves comparison of 
paths from x to y (which exist because the contribution 
to dMSF vanishes if x and y are not connected) . We saw 
above that we can assume that the root points Xi, X2 
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of Ti are connected by E'{TL). Further, there are paths 
from X to y through 'He'(h)i because otherwise 'He'{h) 
is either disconnected from both x and y, or is part of 
a tadpole, and in either case the weight vanishes as we 
saw above. To find the MST path from x to y, the task 
can be broken into subtasks, and one of these is first 
to find the MST path through He'Ch) between its root 
points. If the MST path from x to y passes through 
'HE'iH): the portion within 'He'{h) must be this MST 
path. We now show that (within the sum defining dMSp) 
this path 7mst('^e'(h) !''■£'(«)) between Xi and X2 must 
pass through the most costly edge tu of H.. For suppose 
that tu G E'CH), but the MST path does not traverse 
it. Then there is another edge set which is the same as 
E' except that en is omitted, and these terms cancel in 
pairs (note that the MST paths are the same for these 
edge sets). But the terms with en £ E'{Ti,) and e-j-c on 
ImstO^-e' {n)\''^E' (n)) do not cancel in a similar way, as 
removing e-j-i from this edge set leaves the root vertices xi 
and X2 disconnected, and we know that those edge sets 
cancel among themselves. The reason the root vertices 
become disconnected on removing en (so Ti-E'^n) is not 
IPI) is that if not, then a less-costly path (in the sense of 
the ordering -< in section IB 21) between the roots would 
exist. 

It follows that in comparing possible MST paths on 
Qe', the subgraph Ti.E'{n} can be replaced by a single 
edge from Xi to X2 with cost . We use this result to 
define the induced ordering ttq jn for the quotient graph 
for such a 2-point subgraph H; this ordering gives the 
ordering for any diluted edge set E'(Q/T-L). Note however 
that here we are forced to view H as replaced by an 
edge, not a vertex, in the quotient graph. [Further, Ti, 
is bordered by two other edges, and these three edges 
form a chain, which by the general elementary arguments 
given earlier can be replaced by a single edge of cost the 
maximum of the costs of the three edges, for the purposes 
of finding the MST path 7mst('H_e' kfi') from x to y.] We 
can summarize this whole argument as showing that the 
indicator function can be written as 

4"^[7MST(eB'k£')] = 4""^[7MST(^B'(W)kiJ'(W))] 

X l'f^[lMST{QE'l'HE'iH)WE'(g/n))l (D-27) 

while the edge subsets E'{0/7i) that have to be summed 
over are subsets of the set (£" — E'{H))Li {en} (the latter 
change cause no difficulty, and again the three edges in a 
chain can be replaced by one, with the cost as described 
above). This then shows that the weight factors as 

rfMSF(^|7r) = dMSF{T-C\TrH)dMSF{G /Hlng/n), (D.28) 

where the cZmsf for TC is that for the path from xi to 
X2 to pass though en, while the second simply requires 
a path on the quotient graph GfH to pass through z. 
However, the argument already given above for the MST 
path within Ti on the diluted edge sets E'{7i) shows that 
if the former condition is dropped, then the evaluation of 



the sum is the same. That is 

rfMSF(^kw) - d{n) (D.29) 

for dusFi'HlTTn) with the MST path vertex at the most 
costly edge of H under the ordering nn ■ Thus we obtain 
factorization in the same form as before, as desired (see 
Fig. O. This agrees with the result for an ordering in 
class 7r[o] (because then irg^n = ''^Q-h), but gives the 
correct generalization to other orderings, for the case of 
a 2-point subgraph. For other orderings, the highest cost 
in Ti. has to be compared with those in the remainder 
of G', we emphasize again this aspect of the definition of 

In the present case, the argument can simply be used 
again if GfH contains a 2-point subgraph. For 3-point 
subgraphs, we expect that a more complicated general- 
ization exists, but we have not looked for it. 




FIG. 5: The factorization property of dMSF{G\Tr) proved in 
the text for arbitrary ordering of edge costs. 



5. Proof of renormalizability 

In this section, we assemble the preceding results to de- 
scribe the divergences of the diagrams or subdiagrams in 
MSF(p) theory, and compare them with the correspond- 
ing ones in the perturbation expansion for percolation. 
Here by a diagram, we mean the corresponding Feynman 
integral, including the sum over orderings of the diviSF^Tr 
factors inside the parametric integral. The initial results 
provide the direct motivation for the renormalization of 
the perturbation series. Then we describe the proof of 
renormalizability. 

We saw in section ID 31 that for a subdiagram H and 
an ordering in the class 7r[o] (or for the whole diagram 
G, and any ordering), the superficial degree of divergence 
wmsf(H, t/|7r) is the same as Lo{Ti). We consider only 
2- or 3-point subdiagrams, including the 2-point sub- 
diagram that contains the path vertex. For any fixed 
ordering on the edges E{G) — E{TL) not in 7i, we can 
consider the sum of dMSF{G\T^)FT^{Ag) over all the or- 
derings T^n of edges in E{H) such that the ordering 
of all edges is in 7r[o] . For each ordering in the sum, 
we saw in section ID 41 that the dMSF weight reduces to 
dMSF{'H\TTn)dMSF{G /Ti-lTTg-n) (for the path vertex case) 
or d{'H)dMSF{G /Hlng-n) (for the other cases). For the 
latter cases in which the path vertex is not in Ti, the 
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weight dMSF(0|7r) is independent of the ordering tt^^, and 
the sum over the latter can be performed using eq. (|D.16p . 
which shows that the part of the associated with the 
subdiagram has reduced to unity, as in eq. (jP.lsp . The 
remaining factor on the right hand side of eq. (jD.16p de- 
pends on the a's for the subdiagram, but of course not 
on their ordering. To leading order as all those a's are 
scaled to zero, the resulting subintegration has exactly 
the same divergence (not only degree of divergence) as 
the corresponding subdiagram in the field theory of per- 
colation near criticality [including the d{Ti.) weight], and 
the remaining F^^ factor is that for the quotient graph, 
FTrg-ni^g /n)- (Note that here we disregarded the pos- 
sibility that H itself contains a subdiagram that is diver- 
gent; this will be handled later.) For the path vertex, 
the leading divergence of the subintegral is not identical 
to any in the percolation theory, as the sum of F^^ over 
TTu does not reduce to a factor unity for the subdiagram; 
the 2-point vertex function with a mass (0^) insertion, 
which it resembles, is different, though it has the same 
degree of divergence [an example of such an integral was 
discussed in eq. ()III.20P ]. Nonetheless, the weighted sum 
of Ftt has similar factorization properties. 

Motivated by these observations, we aim to prove that 
our perturbation expansion for the path vertex function 
can be renormalized in a manner very similar to that for 
percolation. Indeed, the mass, field, and coupling renor- 
malizations will be exactly the same as in percolation, 
even when they occur inside a IPI diagram for the path 
vertex function (we saw earlier that this is so for the self 
energy parts outside the IPI path vertex function, that 
is connected to this function by a single line). For the 
path vertex itself, the renormalization works and takes 
a similar form as that for a mass insertion in percola- 
tion, but the coefficients are different. We will prove this 
to all orders in perturbation theory. First, we will es- 
tablish that it is possible to perform subtractions as in 
ordinary field theory Feynman diagrams (e.g. for perco- 
lation), with the result that our subtracted amplitudes 
are non-diverging for each graph in every order in per- 
turbation theory. Then we will show that, because of 
a cancellation of some sub-leading pieces involving self- 
energy insertions, the subtractions take the same form 
as in percolation, as indicated above. This then leads al- 
most immediately to the RG equations, and the epsilon 
expansion for the exponents. 

The idea for rendering the Feynman integral associated 
to a diagram finite is intuitively simple. We identify all 
the diverging sub-integrations associated to subdiagrams 
of the types already listed above (called "renormaliza- 
tion parts"), for which the divergence is related to the 
behavior of the integrand as a corresponding set of a pa- 
rameters is scaled to zero, and then subtract away these 
parts of the integrand. One would hope that the resulting 
integrand is then convergent, and even absolutely con- 
vergent. It is necessary to prove this non-obvious result, 
which we will do using results from the literature. 

The procedure is somewhat complex because a given 



diagram may contain several diverging subdiagrams. 
The subdiagrams may themselves contain diverging sub- 
subdiagram (these are revealed by considering several di- 
latation parameters pi attached to distinct subdiagrams, 
which go to zero in some order). These possibilities cause 
no problems for disjoint subdiagrams (that have no com- 
mon edges), nor for a nested sub-subdiagram (entirely 
contained in a subdiagram). The case of subdiagrams 
that are neither disjoint nor nested, called "overlapping 
divergences", is more difficult, but turns out not to be 
a problem. One makes subtractions corresponding only 
to non-overlapping and nested subdiagrams. The pro- 
cedure was defined by Bogoliubov and Parasiuk in re- 
cursive terms [i^, finiteness was proved by Hepp 
and a non-recursive definition in terms of "forests" was 
given by Zimmerman 50]. Together, this formulation is 
called the BPHZ method. These authors worked in terms 
of momentum-space integrals. For the later formulation 
and proofs within the parametric formulation, sec Refs. 
HHiS and the review in Ref. [3|. 

We briefiy outline the result due to Bergere and Lam 
[s?! that we will use. First it will be useful to introduce 
the "generalized Taylor operators" T" [36|, (S^]- For a 
function f{x) of a positive variable x that behaves as 
f{x) ^ aox" as a; — > (ao 7^ 0), such that x"" f[x) 
is infinitely-diffcrcntiable on [0,a) (a > 0), and for our 
purposes with v an integer, (such a function is said to 
have the Taylor series property) these are defined for any 
integer n to extract the Laurent-like series of terms: 

r"/(2;) = a^x" -f aix'^+i + . . . + a^+„a:", (D.30) 

(where oq, ... a^+n are complex numbers) with properties 
T"/(a;) = if n < i^, and (1 - T'')f{x) ~ a;' with 
q > n. While the series has the Laurent form, we do not 
assume / is complex differentiable away from 0, and the 
coefficients can be calculated from / at positive x only, 
by ordinary Taylor expansion of x~'^ f(x) at a; = 0. For 
a function of several variables xi, a;2, . . . , we may define 
generalized Taylor operators 7^"' similarly by acting with 
one of them at a time, but we must be careful as they do 
not generally commute. 

In the following these operations will be applied acting 
on some subset E' of a's for a graph Qhy a. dilatation pa- 
rameter p as p — > 0, and then setting p = 1 in the result: 
TS,f{Ag) = V/({p2a. : i G E'},{a. : i ^ E'})\p^,. 
Thus these extract precisely the leading and subleading 
terms that we have been discussing, up to order n. Here 
when /({p^Q!i : i € E'},{ai : i ^ E'}) has the Taylor 
series property as a function of p, we say it has it with 
respect to the set E', and it is in this case that the oper- 
ator T^, is defined. 

We will need some definitions for properties of the func- 
tions to which the Theorem applies. We will consider 
what Bergere and Lam [s^] call a "nest" of edge subsets, 
which is a filtration, that is a, set Af — {Ei, . . . , Er} of 
edge subsets such that 

C Si C £^2 C • • • C £:* C E{g) (DM) 
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in which the inclusions are strict except possibly the 
last. For a function Z{Ag), we say that it has the "si- 
multaneous Taylor series property" with respect to the 
filtration M if there is a set of integers i^e,. such that 
(n*=i Pr ) Z {p'^ Ag) has simultaneous Taylor series in 
the set of pr near, and does not vanish at, Pr — for all 
r; here p^Ag stands for the ordered set Ag of a's, but 
each Ui acquires a factor p^. for each subset Er to which 
i belongs. For example, the function l/(ai + 0^2) has 
the simultaneous Taylor series property for the filtration 
i?i = {l}, i?2 = {l,2}. 

Now we can state (a special case of) the theorem of 
Bergere and Lam: if (i) Z{Ag) is infinitely differentiable 
for < ai < 00; (ii) Z{Ag) and its a derivatives are 
polynomially bounded when arbitrary subsets of Ag are 
scaled to 00; (iii) Z{Ag) has the simultaneous Taylor 
series property with respect to every filtration M of edge 
subsets, then the integral 

Ir= it da,e-^-"'*i?[Z(Ae)] (D.32) 

with i > is absolutely convergent. Here the R operation 
is the subtraction operator which can be defined as 

AA E'eAf 

where the sum is over all filtrations JV of the set Ag of 
a's. 

In its general form, the theorem applies to many inte- 
grals that are not related to Feynman diagrams in any 
obvious way. Now we wish to apply it to the Feynman 
integrals in our perturbation expansion. First we point 
out that these integrals do satisfy the hypotheses of the 
theorem. Indeed, the integrands of our integrals contain 
factors that occur in the field theory of percolation, which 
for this purpose is no different from a cubic-interaction 
scalar field theory, times the factor F,^ for some order- 
ing TT (times c?MSF and summed over tt, but we need not 
consider this here; this sum can be exchanged with the 
integral and then taken under the R operation if desired) . 
The integrand in the cubic theory satisfies the conditions, 
and it is easy to see that the i^^ factor does not change 
this. 

To go further, we note that when applied to Feynman 
integrals based on a graph (the graph made no appear- 
ance in the statement of the theorem) , R can also be ex- 
pressed in many other ways, one of which is as the sum 
over forests of renormalization parts [37| . As we know, a 
forest is a collection of trees, but here the trees are not 
spanning trees on our lattice or our graph Q. Instead, a 
forest is any set of renormalization parts in Q (which are 
IPI connected 2- or 3-point subdiagrams), such that for 
any two such parts in the set, either one is entirely con- 
tained in the other (both for its vertices and its edges), 
or else they are disjoint. (Often in the literature, a for- 
est is pictured as a set of non-intersecting boxes overlaid 



on the depiction of the Feynman diagram.) In this form 
for R, the sum over all filtrations is replaced by a sum 
over all forests, and each edge set E' in the product is 
that of a single renormalization part belonging to that 
forest. We note that in the BPHZ formulation, whether 
in parametric form or not, no divergent integral or cutoff 
is mentioned. The subtractions are carried out instead 
on the integrand (which however, before the subtractions 
are performed, does have the property of diverging more 
strongly in some limits). 

We will apply the Theorem to the integral for a di- 
agram Q, in which the sum over orderings, and dMSF 
factors, are taken into the integrand. That is, 

Z{Ag) = 

TT 

X p-^/^Ag)eM-'i^LA'''^e.t)- (D.34) 

Our earlier remarks imply that for each renormalization 
part, the subtractions (in forest form) remove precisely 
all the superficially divergent pieces and no more. Thus 
in this form, the R operation is exactly what one might 
expect it to be from the discussion preceding the state- 
ment of the theorem, and the theorem says that these 
subtractions result in an absolutely convergent integral. 
For subdiagrams H of a diagram Q, these subtractions are 
exactly the same as those for a diagram of the same type 
(number of external points, and presence or absence of 
the path vertex), with one exception. This is the sublead- 
ing superficial divergence in the case of 7i a self-energy 
subdiagram. In the subleading generalized Taylor expan- 
sion (that is, r-2|£^(«)l _ 7— 2-2|£;(w)| acting on the di- 
latation factor for the subdiagram), part of it comes from 
expanding F,^ to order a (for some a in the subdiagram) 
times the leading term from the percolation integrand; 
the terms from F^^ arise from orderings in class tt^j^j and 
from subleading terms in class ttjq] . This does not cor- 
respond to the subtraction made to any whole diagram, 
and would thus be difficult to include in the renormaliza- 
tion scheme. Fortunately, these subtractions cancel, not 
for the given diagram, but between diagrams of the same 
order that differ only in the placement of the self-energy 
insertion in the graph. This cancellation result will be 
proved in section fP 61 below. 

Hence because we are always interested in the sum of 
all diagrams in each order anyway, the only subtractions 
that have to be made correspond to those that would 
be made to Q when it is a renormalization part. It fol- 
lows that the subtractions correspond to subtracting the 
Taylor series in for the subdiagram, where k is the 
wavevector entering the subdiagram, and replacing the 
original graph by the quotient by the subdiagram, times 
these Taylor coefficients in place of the subdiagram. For 
the vertex and path vertex cases, the subtraction is sim- 
ply at zero wavevector, while for the self-energy the first 
order term in k^ has to be subtracted also. This is easily 
seen, as the subtractions to Q itself are just its Taylor 
expansion in k^ to the given order [H, [H, [H, [s^]. If 
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we include the zero-loop parts of the vertex functions, 
this implies that the renormalized vertex functions in this 
renormalization scheme obey 



.(2) 



k=0 



dk2 « 

rf {0,0,0; g,t) 
rg''^(0,0,0;g,i) 



rg'P^)(0,0,0;g,i) 



t, 
1, 

1, 
1. 



(D.35) 



In view of the condition on the 2- and 3-point vertex func- 
tions at k = 0, the coupling and mass-squared appearing 
in the propagators in the expansion can be identified with 
the renormalized values, so there is no subscript zero on 
these quantities. 

Now that the renormalized perturbation series defining 
the Tji are known to be finite, for example at non-zero 
wavevectors away from the point at which the above 
conditions are given, we can modify the renormalization 
scheme. Namely, we can add a finite part (more accu- 
rately, a series of finite terms) to each subtracted piece 
in the definition of the renormalized integrand. These 
can be chosen in each order to modify the renormaliza- 
tion conditions, and the combinatorics again works out. 
This changes the renormalization scheme, and for exam- 
ple we can modify the conditions above to specify values 
at non-zero wavevectors (except for the mass-squared): 



= 1, 



|q| — K 



rf{{q.}; g,t) 


SP 


= 9, 


(D.36) 
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= 1, 






SP 


= 1. 





Here SP — symmetry point denotes a symmetric con- 
figuration of external momenta qi, q2, qs, which (by 
rotational symmetry) we take to be any triple satisfy- 
ing q2 = ^2 for i = 1, 2, 3, q, • q^ = {i ^ j). 
Note that g and t now have a different meaning than 
before. In this form, we can now set t — and work 
directly at the critical point, as in each order in pertur- 
bation theory the non-zero wavevector scale k prevents 
the left-hand-sides from diverging in the infrared (the 
self-energy ~T^^^ {0; g,0) is not infrared divergent). This 
renormalization at zero mass-squared is quite convenient 
technically. 

As we mentioned above, the BPHZ subtraction scheme 
requires no reference to, nor use of, a cutoff. It is possible 
to develop the RG equations directly from this scheme, 
working with non-divergent expressions only, and lead- 
ing for example to the Callan-Symanzik equation when 



the renormalization scheme at zero wavevector, non-zero 
t is used However, for calculational purposes, we 

prefer to write intermediate quantities in terms of expres- 
sions that diverge as A — > oo as in traditional approaches. 
The bare vertex functions are given by the original, un- 
subtracted Feynman integrals with cut-off, including as 
always the rfMSF-FV factors. For emphasis, we now write 
these as Fq's. They are viewed as functions of the bare 
coupling go and mass-squared to ; as well as the wavevec- 
tors and cutoff A. Then all the subtractions that de- 
fine the renormalized amplitudes can be collected into 
changes of the parameters to g and t — 0, and changes 
in the scale of the "operators" (p, (fP' and that described 
by the path vertex. That is 



Z^'\go, «, A)Z^. (go, «, A)F[,"^"^({qa, {q,}; 50, to, A) 



{NX) 



= F 



(N,L] 



({qJ,{qi};5,«^,A), 



(D.37) 



Z0(go,'«,A)Zpv(go,t,A)F^ ' ^{qJ; ffo, to. A) 

= F(,''P^)({qa;g,^,A), (D.38) 

and in the limit A — ^ 00 the dependence of all F/? on 
A drops out. These equations require five equations to 
define the dependence of Z^, Z^2, Zpy, go, and to on 
g, K, and A, and these are provided by the five condi- 
tions (|D.36p . when these are expanded in perturbation 
theory in go. At this point the treatment of our theory 
has come to closely resemble an ordinary field theory, the 
main difference being the form of the Feynman rules for 
calculating F'^^.pv)^ ^pj^g j-^ggi; important conclusion of 
the analysis is that the path vertex is renormalized mul- 
tiplicatively by Zpy. We describe in the main text the 
derivation of the RG equations, and the calculation of 
exponents to one-loop order. 



6. Cancellation proof for subleading terms 

In this section we present the proof that the particular 
subleading terms in the Laurent expansion as the a's in 
a self-energy part (not containing z) go to zero, that do 
not appear for the self-energy in an external line, actually 
all cancel among graphs with the same self-energy part 
inserted in different edges. 

First wc show that the dependence on the a's in a self- 
energy of the weighted sum of factors simplifies. We 
suppose throughout this section that we consider a fixed 
graph Oo with an ordering ttq, and we then modify this 
graph to obtain Qi by inserting a given self-energy graph 
H on an edge i of Go- Thus, the edge i is replaced by 
two edges i', i" , with the self-energy 7i in between. In 
the parametric integral for the diagram, the parameter 
ai in Qo is replaced by parameters a^' , a^" , and there are 
additional parameters for the edges of TC. 

We know from Section ID 41 that the diviSF weight for 
Qi is determined by the ordering TTg-fu in which the cost 
replacing the original £i is the largest of ii' , £i>' and those 
in Ti, independent of how these are ordered relative to 
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each other. Moreover, the weight factors as 

dusFiQilT^) = dMSF{'H)dMSF{Gi/'H\TTg^/T^). (D.39) 

Throughout the argument, we wiU compare cases in 
which TTg-f-H = ttq (in an obvious sense), and is fixed. The 
sum of the F-^ factors over the orderings of the edges that 
replace i can be calculated, and this is done most easily 
by returning to the original calculation of F,^ from the 
action of the Omsf operator in section ID^ The desired 
sum has the effect of simplifying the integro-differential 
operator to the following form, and acting under the 
parametric integral gives 




(D.40) 



where (i) D is the usual (|i?(5o)|-dimensional) integration 
domain for Gq with ordering ttq, defined by > • • • > 

i7ro(|-B(c?o)l '■0' the integrand, t^/, ti" , and the i's 

associated to edges in H are all set equal to ti, and (ii) 
are the same as aj except for a^, which is the sum of a^' , 
ai" and all the a's in Ti.. The product on the right-hand 
side is simply F^^ for Qq but with this substitution; we 
denote it Fl^^{i). 

A second trick that is commonly used for parametric 
integrals is also useful: if the integrand only depends 
on the sum of two parameters, say ai and a2, then these 
integrations can be combined into a single integral over a 
which takes the place of ai+a2, at the cost of introducing 
a factor a into the integrand: 

J daida^ . . . — j daa . . . . (D.41) 

(This can be generalized to any number of variables, but 
we do not require that.) It can be shown that, like the 
weighted sum of factors, the rest of the parametric 
integrand only depends on the sum ai' + ai" (this can 
be shown by some further use of the relation, mentioned 
earlier, of this integrand to combinatorics of weighted 
spanning trees, which we do not enter into). Then we 
use this result to recover an integral over a single ai in 
place of those two. 

We now consider the generalized Taylor expansion of 
the integrand with respect to a dilatation parameter p 
applied to the a's in H. The leading behavior is seen to 
give simply the integrand for Ti. from percolation (with 
zero external wavevectors) , times the integrand and -FV^ 
factor for the quotient graph, as discussed above. We 
now turn to the subleading terms of a particular form: 
those that come from expanding the above factor F^^ (i) 
to first order in (or simply in a's in Ti), times the lead- 
ing behavior of the rest of the integrand, as p ^ 0. The 
rest of the integrand factors into that for the quotient 



graph times that for the subgraph, and we make use of 
the technique for replacing the a^/, by a^. This factor 
in the integrand is now independent of which edge i of 
Qq was chosen for insertion of Ti. The dMSF factor is also 
independent of i because of the choices of ordering made 
earlier. Let us write ^ as.c. for the sum of a's associated 
to Ti. There are two types of terms in the expansion 
of F'^^{i) at first order in X) Q^s.c. for each i: a) those in 
which ^ as.c. appears in the numerator, in which case it 
replaces a^, so giving the factor 'Y^a^.e./ ai times b) 
those in which it comes from expanding a denominator, 
which must be one of those indexed j > nQ^li). This 

gives a factor - J2 cts.e./ J2k=i Q^jroCfe) times F^^. We re- 
member to multiply by ai (because we replaced two a's 
by this one), and then sum over the positions i of the 
self-energy insertion. This gives 

= 0, 

times other common factors. This completes the proof. 

What we have shown is that this type of subleading 
term actually cancels in the sum of diagrams in each or- 
der. Because the superficial degree of divergence of a self- 
energy diagram at six dimensions is 2, these subleading 
parts are also superficially divergent, and are subtracted 
for each diagram by the R operation defined earlier. We 
view the present result as showing that the subtracted 
terms cancel, and because the sum of the subtracted in- 
tegrals is finite, these subtracted terms can be dropped 
and the result is still finite. The cancellation is indepen- 
dent of other subtractions associated with renormaliza- 
tion parts, which might appear either inside Ti or disjoint 
from it. In particular, this allows the cancellation to be 
made for any number of self-energy insertions in a graph. 
The remaining terms subtracted from a self-energy inser- 
tion are then exactly those that occur in percolation, for 
any number of self-energy insertions. 

The principle underlying these pleasant cancelations is 
not entirely clear to us. It seems likely (because they in- 
volve the derivative of the self-energy in percolation with 
respect to t) that they are associated with the notion 
that the Omsf operator should be renormalized, so that 
it acts on t rather than to- This operator has the prop- 
erty that it is invariant under any reparameterization of 
the variables ^ — T{t^) provided that T is a mono- 
tonic function and has no explicit dependence on e (this 
invariance of the geometry of MSTs is related to that 
emphasized in Ref. 12\). This property of Omsf was not 
explicitly used anywhere in our construction. We will 
not attempt to give here a conceptual proof using these 
ideas. 
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